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Phase transition is a war game. It widely exists in different kinds of complex system beyond 
physics. Where there is revolution, there is phase transition. The renormalization group transfor- 
mation, which was proved to be a powerful tool to study the critical phenomena, is actually a game 
process. The phase boundary between the old phase and new phase is the outcome of many rounds 
of negotiation between the old force and new force. The order of phase transition is determined by 
the cutoff of renormalization group transformation. This definition unified Ehrenfest's definition of 
phase transition in thermodynamic physics. If the strategy manifold has nontrivial topology, the 
topological relation would put a constrain on the surviving strategies, the transition occurred under 
this constrain may be called a topological one. If the strategy manifold is open and noncompact, 
phase transition is simply a game process, there is no table for topology. An universal phase coexis- 
tence equation is found, it sits at the Nash equilibrium point. Inspired by the fractal space structure 
demonstrated by renormalization group theory, a conjecture is proposed that the universal scaling 
law of a general phase transition in a complex system comes from the coexistence equation around 
Nash equilibrium point. Game theory also provide us new understanding to pairing mechanism and 
entanglement in many body physics. 
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1. INTRODUCTION 

When physicists encounter millions of interacting 
molecules or atoms, they can not control the trajec- 
tory or momentum of each individual particles, so they 
study the macroscopic states at different temperature or 
other physical parameter. In most cases, there would be 
some significant change of the macroscopic states at cer- 
tain value of parameters, they call it a phase transition. 
Phase transition are common phenomena in all branches 
of physics. People's interest in phase transition can trace 
back to thousands of years ago. A recent example is the 
superfluid-Mott-insulator phase transition occurred in a 
gas of ultracold atoms in optical lattice [![. 

Statistic mechanics were developed to provide a theo- 
retical description of phase transition. The occurrence of 
phase transition is related to the singularity of statisti- 
cal functions in the thermodynamic limit (see Ref. Q for 
review). However statistical mechanics is not powerful 
enough to predict all different kind of phase transition. 
For classical Hamiltonian systems, the hypothesis con- 
necting phase transitions to the change of configuration 
space topology was proposed recent years (see Ref. [3( and 
references therein for review). Topological order also rise 
in quantum systems (see Ref. |J| for review), such as 
the fractional quantum Hall system. Different fractional 
quantum Hall effect states all have the same symmetry, 
that is beyond the Landaus symmetry breaking theory. 
The topological order in quantum system is related to 
degenerate ground states, quasiparticle statistics, edge 
states, or momentum topology [a], et al. These topolog- 
ical order theory told us many interesting phenomena 
in some special models and special systems. However 
it is still far from providing us an universal explana- 
tion to phase transition occurred in all different branches 
of physics, much less in other complex system beyond 
physics. 

Generally people believed that it is the thermal fluctu- 
ation that drives the transition from one phase to another 
in classical phase transition. As temperature is lowered, 
the thermal fluctuations are suppressed. The quantum 
fluctuation began to play a vital role in quantum phase 
transition. Unfortunately this kind of argument does not 



hold in many quantum system. More over, some physi- 
cist believed that phase transition is due to the competi- 
tion between two competing orders in a physical system, 
but people can always propose some anomalous examples 
which can not be explained by two competing orders. 

The occurrence mechanism of phase transitions is still 
not clear in countless systems. Similar critical phenom- 
ena arose in a broad physical and social systems. Many 
unrelated models covering physics, chemistry, biology 
present similar scaling laws. Anyone who saw this can 
not help asking why. Is there a 'theory of everything' 
that can give us an universal explanation? Such a theory 
sounds like superstring theory. 

What I am trying to do in this paper is not to establish 
the superstring theory of phase transition, but to present 
an universal theoretical explanation to phase transition 
occurred in different systems based on renormalization 
group theory and game theory and beyond the two, such 
as topology, quantum field theory. 

The first step is to break the envelop of physics and ex- 
tend the concept of phase transition from physical system 
to complex system including chemistry system, biologi- 
cal system, social system, economic system, et al. Phase 
transition is a sudden jumping from one stable state to 
another. A two player game always has two stable states, 
the winning state and losing state. Thus we can define 
phase transition as a war game. If we check the war game 
carefully, one would see it has all the the same phenom- 
ena as phase transition occurred in physical system. War 
is a conflict between millions of soldiers who are armed 
and well organized. The butterfly effect is a basic char- 
acter of war going on. The final destiny of the fighter is 
determined by some minor accidental event. When the 
two groups with equivalent force are fighting against each 
other but keeping at a draw state, if any one of them is 
reinforced a little, he will win the war in a few seconds. 
This is a phase transition. 

So we can take phase transition as a war game, the 
strategies of the players extended the strategy base man- 
ifold. When we are studying the state evolution of the 
game corresponding to different strategy the topology of 
the strategy base manifold comes in. The phase tran- 
sition is a war game between new phase and old phase, 
each of them is governed by a kind of dominant inter- 
action^ may has many minor affiliated companions). It 
will be shown in the main tex of this paper, the surviving 
strategy of the two phases carry opposite winding num- 
ber, the sum of these winding numbers is a topological 
number on the strategy base manifold. After the wind- 
ing number are annihilated by pairs, the last one winding 
number around the last surviving strategy is decided by 
the topological number, this also decided who will be the 
winner between the old phase and new phase. It is in 
this sense, we call it a topological phase transition. In 
fact, phase transition is always related to the topology of 
the strategy base manifold when the strategy base mani- 
fold is compact. In some cases, the base manifold is open 
and noncompact , there is no need to consider the effect 
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of topological constrain. In fact, the compact strategy 
manifold like a finite region confined by fencing, it put 
some constrains on the way we choose strategy. There is 
singular point we can never eliminated smoothly in strat- 
egy space, that is the fundamental origin of topological 
phase transition. 

The paper is organized as follows: 

In section 2, the most general conception of phase tran- 
sition in complex system is defined. 

In section 3, we established the game theory of renor- 
malization group transformation, and find the general 
solution of renormalization group transformation equa- 
tions. The fundament classification of phase transition 
through symmetry losing is presented. 

In section 4, topological current theory of phase tran- 
sition is established, this theory spontaneously produced 
an universal equation of phase coexistence. A conjecture 
on the universal scaling law is proposed base on a topo- 
logical hypothesis in fractal strategy space of game the- 
ory. Further more, we established the evolution equation 
of phases and the quantum phase coexistence equation. 

In section 5, we developed the quantum statistics of 
many player game, and proposed a conjecture to find 
the fixed point of a many player game using quantum 
density matrix. More over, a new quantity to measure 
the entanglement of quantum states in a game is found. 
The single direction of renormalization group follows the 
second thermodynamic law, the renormalization flows to 
increase the entropy as well as the quantum entangle- 
ment. More over, the coexistence state of multi-player 
game is discussed. 

In section 6, we developed the game theory of phase 
transition in classical many body system as well as quan- 
tum many body system. A new holographic topolog- 
ical quantity to characterize momentum space is pro- 
posed. We studied the quantum many body theory of 
war game and gave a new pairing mechanism base on 
prisoner dilemma. 

The last section is devoted to a brief summary and 
outlook. 



2. PHASE TRANSITION AND WAR GAME 
2.1. Phases of complex system 

A complex system consists of many different elements 
that are connected or related, it appears like a black box 
to us. One can obtain the information within the complex 
system by its responses to external perturbation. These 
responses and perturbations are macroscopic variables. 
Different stable phases of complex system are character- 
ized and distinguished by these observablcs. 

For most condensed matter physicists, a liter of wa- 
ter is a complex system, since it is hardly possible to 
find exact analytical solution for the motion of 10 27 H2O 
molecules. We can characterize its different phases by ob- 
serving its chemical composition and physical properties, 



such as volume, pressure, temperature, density, crystal 
structure, index of refraction, chemical potential, and so 
forth. One of the main task of physicist is to find the rela- 
tions among these macroscopic variables via presumably 
well known microscopic interactions between particles. 

A stable phase has self-restoring ability when it is dis- 
turbed from equilibrium states. During a relative long 
life time, the chemical makeup does not decompose, and 
the physical properties keep a good manner. If we know 
the position and momentum of every particles in a dy- 
namic system at a given time, the evolution of an integral 
system is exactly predictable. Unfortunately, there is few 
integrable many body system. The position and momen- 
tum of particles can not be exactly measured at the same 
time due to the Heisenberg uncertainty principal. So we 
define the state vector of a stable phase by the complete 
set of observables. 

Let x be n independent states of a complex system. 
For a physical system consists of N interacting particles, 
the components of an arbitrary vector Xi are consists 
of 2>N position coordinates <?i, <?i, <?3jv, 3iV momentum 
coordinates ...,P3N, and a vector of other physical 

parameters 7, such as temperature 71 = T, pressure 72 = 
P, particle density 73 = N, volume 74 = V, chemical 
potential 75 = fi,- ■ ■ , conductivity 7^ = a, susceptibility 
7fe = Xi an d so on. Here the state vector x = (q,p) is 
a much more general conception than that of statistical 
physics. It indicates the information inside the black box 
of any complex system. 

The response of the complex system are induced by 
external input vector 7. For physicist, 7$ represents 
those familiar external applied magnetic field, electric 
field, pressure, neutrino current, electric current, detect- 
ing laser beams, heaters, and so forth. For chemist or 
biologist, the input vector represents something like en- 
zymes, chemical accelerator. 

Not all the components of the state vector are observ- 
able under perturbation. The output vector y represents 
those observables that people can definitely detected in 
laboratory. The output is strongly depend on what peo- 
ple want to study. For example, there are circulatory 
systems, nervous systems and digestive systems within 
human body. If we are studying the human population, 
there is no need to take into account of these subsystems; 
one only counts the people, the output vector covers the 
number of people, distribution of people, etc. If the sub- 
ject is about flu's spread, it may be best to discuss the 
immune subsystem. 

For a general dynamic system 

= F in (x,j), O ut = F out (x,j), (2.1) 

the output vector is not always diffcrcntiable for all de- 
grees of differentiation on the whole range of the param- 
eter space M(a?, 7). There exist critical points (of, 7*) 
at which the C k output functions blow up. The whole 
input space is divided into separated blocks by these crit- 
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FIG. 1: The whole phase diagram is split into discrete do- 
mains. The stable phase exist in the inner region. The do- 
mains wall is the phase coexistence region. In each domain, 
there is a dominant interaction which is the king ruling over 
the other weaker interaction. There is a war going on at 
the phase coexistence boundary. When we tune the physi- 
cal parameters, we are helping a certain interaction to fight 
against the others, this certain phase would grow stronger and 
stronger, it will finally unified the whole phase diagram. 



ical points, 

M°(r) = {(rJ.rS) U (r*,r*) • • • U «,<)}, (2.2) 

here we denote (r :— (x,j)). The stable phase are defined 
in these discrete blocks. The output functions present 
very good behavior within the blocks but diverge at the 
very boundary. These boundaries are where the phase 
transition occurs. 

This general mathematical definition of phases for 
complex system applies for many different fields. The 
most familiar Ehrcnfcst's dcfinition[6( of phase transition 
in thermodynamics is a good exemplar. The output vec- 
tor is only a function — free energy. The input vectors are 
thermodynamic quantities, temperature T and pressures 
P. For the zeroth order stable phase the free energy of 
the two phases F(T, P) is C°° in the whole input space 
(T,P) E [— oo,+oo]. For the first order phase transi- 
tion, the free energy F(T, P) is continuous in the region 
T G [T h T f ], P G [P,,Pf], but the first order derivative is 
not continuous, 

dF A dF B 8F A dF B 

g T T QT > gp r Qp > 

T A G (Ti,T c ),T B G (T c ,7», 

P A &(Pi,Pc),PB&(Pc,Pf)- (2.3) 

The first order stable phase blocks are divided into 
smaller blocks by the second order phase transition. 

The phase transition of complex system may be defined 
from the divergence of the C k output functions. They 
could be any observable functions. For mathematician, a 
stable phase is marked by a C k function in the discrete 
blocks on the input vector space. The critical point is 
the phase boundary between the separated blocks. 

2.2. Phase transition 

The stable output states of a dynamic complex sys- 
tem are confined in different domains in the whole input 



state vector space. As the input vector changes within 
a domain, the system is in a stable state, it represents 
a kind of physical order. When the input vector jumps 
from one domain to another, the system jumps from one 
stable phase to another. The domain wall is the exact 
phase boundary at which the old phase becomes unsta- 
ble and decays, but the new phase comes into being and 
finally leads to new stable order. Phase transition occurs 
everywhere in nature. Every phase transition indicates 
a revolution induced by the interaction between the sys- 
tems and their environment, the fittest states of the old 
domain is replaced by the fittest states of the new do- 
main. 

The lifeless nature world is much more intellectual 
than most physicists thought. What physicist measures 
in experiment is always the observable in equilibrium. 
Quantum theory tell us the energy levels of molecules 
and atoms are discrete and quantized, we can calculate 
the transitions probability between those levels involv- 
ing the absorption or emission of photons. But we don't 
know how and why that occurs. In principal, the non- 
equilibrium process of many particle system follows the 
same rules that govern the living systems, such as ants 
society, honeybees group, traffic system, etc. 

If we focus on the behavior of particles before they 
reach equilibrium, one would doubt that electrons are 
possibly intelligent particles. One simple example is two 
resistors connected in parallel in electric circuit. As all 
knew, the current is proportional to the inverse of resis- 
tivity in subways. In the beginning, the electron moves 
together in the main path. When they reach the bifur- 
cation point, they split into two subways. Less current 
in the strong resistive way, and larger current in the less 
resistive one. Like the cars in traffic, if all the electrons 
choose the less resistive way, they block each other until 
it becomes more resistive than the other subway. Then 
some electrons withdraw and transferred to the other 
way. There is no traffic jam in electric current network, 
because the cars are everywhere under the control of a 
global electric field. 

The non-equilibrium dynamic process shows up in the 
vicinity of critical point, at which a small quantitative 
change of a parameter would results in a qualitative 
change of the global behavior of a complex system. Bi- 
ological systems that adapt to their environment thrive, 
those that fail to evolve fade away. The environment are 
external input of the biological system. A stable phase 
of a biological system only exist in a finite region of en- 
vironment parameters, so does a lifeless physical system. 
The stable phase of a collection of H2O molecules is wa- 
ter between 0°C~I00°C, below 0°C(273 K) is ice crystal, 
and it becomes gas above 100°C. 

Like any living creatures in nature, a lifeless physical 
system of interacting H2O molecules evolutes following 
the rule — "survival of the fittest". They took differ- 
ent collective structures according to different pressure, 
temperature, impurity, radiation, gravitation field, elec- 
tric field, magnetic field, density, etc. The stablest phase 
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at certain range of the parameter space survives as the 
fittest structure, the other phases in this domain fade 
away. Out of this special domain, the molecules have to 
reorganize their collective motion pattern in order to fit 
the new environment. When the new phase is born, the 
old phase dies. 

The phase transition most scientists observed in 
physics, chemistry, biology, complex network, economics, 
• ■ • is a revolution. Phase transition is a war between the 
force of old phase and the force of new phase. The funda- 
mental particles are millions of armed soldiers. They are 
divided into different large-scale armed groups who are 
fighting with one another. There is always a critical re- 
gion in which the fighters bet their bottom dollar on one 
war, the winners take everything, the losers get nothing. 
This critical point is the phase boundary. Sometimes, 
the turning point is not so definite that it broadens into 
a finite region, this indicates a crossover transition. 



3. RENORMALIZATION GROUP THEORY 
AND GAME THEORY 

3.1. Game theory and Renormalization group 

Wilson's renormalization group theory provides a fun- 
damental non-perturbative approach to quantum critical 
phenomena [7[. One basic character of quantum many 
body system is that the microscopic particles are not 
distinguishable. Army is the best social system for sim- 
ulating collective behaviors of quantum many particles 
in physics. Usually the soldiers are identical particles in 
the eyes of a general, but are distinguishable particles for 
a sergeant. The hierarchy structure of army plays the 
same role as quantum numbers in physics. The statis- 
tics of particles is scale dependent. One example is two 
column of dipolar Bosonic atoms obeying anyonic statis- 
tics due to the long range interaction. There are other 
biology system which similar to human society, such as 
ant group, honeybee group, we mainly take the army as 
a basic example to demonstrate renormalizattion group 
theory. 

In Kadanoff construction^, a certain number of neigh- 
boring particles are grouped into one cell which act as 
new elementary particles of the renormalized Hamilto- 
nian. At critical point, the Hamiltonian is identical to the 
original Hamiltonian. For an army, this coarse-graining 
procedure naturally take place. An army has a hierar- 
chical structure, the units of different size include a col- 
lection of lower rank of subordination particles. 100 men 
are group into a company, each company acts like one 
particle at higher rank marked by a captain. Every 10 
companies form a regiment, the regime particle may be 
named after by colonel. This coarse-graining procedure 
take finite steps from brigade to division, to corps, and 
finally to an army. 

In fact, the coarse-graining procedure is a simplifica- 
tion of army at war. In the microscopic level, it is the 



soldiers who are fighting with each other. Since they do 
not hate each other personally, they behave as indistin- 
guishable identical particles, and fight as a whole. Renor- 
malization theory simplifies the war between millions of 
soldiers to a war between thousands of companies, to a 
war between hundreds of regiments, finally to the war 
between two army. It is the war between two generals, 
it is also the war between hundred of colonels as well as 
captains. 

The mean field theory view the war as a fight between 
two full generals dressed up by millions of soldiers. This is 
correct in most cases, but it is not always accurate, for the 
general's strategy is carried out by hundreds of colonels 
and captains instead of the elementary soldiers. At the 
critical point, the correlation between the members of 
an army extends to its maximal value. As the butterfly 
effect says, the battle may be lost due to a nail which fail 
to fix the shoes of the horse. For want of the horse, a 
rider is lost. The lost of one rider may directly leads to 
the loss of a battle. 

The renormalization group starts from the most fun- 
damental particles of the army: soldiers. The first order 
renormalization is to reduce the interaction between mil- 
lions of soldiers to hundred of captains which are dressed 
up by soldiers. The second order renormalization pro- 
cedure is to identify effectively the captains belonging 
to one regiment with one colonel. This particle-blocking 
process may continue, and finally end up with full gen- 
erals. 

This explanation of renormalization group theory from 
war between armies is not merely a parable. There is a 
rigorous mathematical correspondence between the game 
theory of war and renormalization group theory. Let's 
take the two-dimensional Ising model on triangle lattice 
as one example. The soldiers are spins cr, , the battle field 
is the triangular lattice, the Hamiltonian of this game 
reads 

JV JV 

H = 7i ^ o% u i + 7a ^ gj, ((T 4 = ±1). (3.1) 

<ij> i 

7i = J/ksT denotes the interaction between the nearest 
neighbor spins. 72 — fiB/ksT is the effective external 
applied magnetic field. This model may be treated as N- 
players game, the N spins are N players, each of them 
has two strategy ±1. The Hamiltonian is the payoff func- 
tion. The players take different strategy to minimize the 
energy function. Another different modelling of this Ising 
model by game theory is to take it as a two-player game. 
We christen the two players 71 and 72. 71 's task is to 
choose a strategy I71) in its strategy space {71} to con- 
trol the interaction between neighboring soldiers, so that 
they act following his orders. 71 governs the soldiers by 
either ferromagnetic interaction or anti-ferromagnetic in- 
teraction. 72 commands the spin soldiers to keep strictly 
in the direction of external magnetic field. 71 and 72 
choose different strategies to win this game. 

A decision rule for 71 is a operator / 7l : {72} =>■ {71}, it 
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associates each strategy |72)g{72} of 72 with the strate- 
gies |7i)g/ 71 I72), which may be played by 71 when he 
knows that 72 is playing I72). Similarly, the decision rule 
for I72) is a map /| 72 ) from {71} to {72}- When a pair of 
strategies (|7i), I72)) satisfies 



l7i}e/ 7l |72>, l72>e/ 72 |7i> 



(3.2) 



they form a consistent pair of strategies. The set of con- 
sistent pair may be empty or very large or it may reduce 
to a small number of bi-strategies. The problem of find- 
ing consistent strategy pairs is so-called fixed-point prob- 
lem. We may construct a consistent map of the strategy 
pair, 

V(|7i), |72» G {71} x {72}, /(|7i), |72» := / 7l x / 72 , 

(3.3) 

such that 

/ab|7i> - / 7l / 72 |7i) = l7i) € {71}, (3.4) 
!ba\12) = / 72 / 7l |72) = IY2) e {72}. (3.5) 
Then a consistent pair is explicitly written in the form, 



l7i) 
I72) 



/- 



72 



71 




l7i) 
I72) 



(3.6) 



This decision rule matrix of the two players is just the 
renormalization group transformation matrix 



U = 



U 

U 



(3.7) 



which can be deduced from the decimation procedure of 
renormalization group transformation. The players may 
play many steps to reach an agremcnt. Each time we 
group the spins in a sum by Kadanoff blocks, the origi- 
nal degree of freedom is decimated into the fewer degree 
of freedom. Every block is now a new giant elementary 
particle whose spin is determined by the majority rule. 
The Kadanoff-blocking process is actually a game pro- 
cess. After the first step of Kadanoff-blocking, the two 
player accomplished the first round of game, and their 
possible strategy space is reduced to a smaller one due 
to the information they get from the first round game. 

The rescalcd Hamiltonian is of the same form as the 
original one, 



N'/2 N' 
<ij> i 



(3.8 



The new parameters represent the new strategy of the 
two players for the second round of the game, they follow 
the renormalization group transformation, 



ii = 71(71,72), 72 = 72(71,72)- 



(3.9) 



We denote the strategy space of the two player as K 
{|7i), 1 72)}, the strategy vector of the two players is 



17) 



l7i) 
I72) 



(3.10) 



The game operator is a map from the strategy space to 
itself, 



U: \-f)eK-> | 7 ') e K. 



(3.11) 



This game operator is equivalent to the renormalization 
group transformation. Applying Brouwcr's fixed point 
theorem, we know if the the strategy space is convex 
compact subsets of finite dimensional vector space, there 
is at least one pair of consistent strategy which satisfies 



|7) - U\i). 



(3.12) 



1 7) is the brouwer fixed point, or Nash equilibrium point. 
This fixed is the saddle point of physical observables on 
the manifold expanded by (71 , 72 ) • 

This game operator for the two dimensional Ising 
model may be derived from the recursion relation the 
semigroup transformation, 



7 
f 



(7i,7 2 ) 



7 = (71,72) 



(3.13) 
(3.14) 



here L is rescaling factor of Kadanoff-blocking. We calcu- 
late a statistical physical observable, such as free energy 



F = 



InZ, Z = Tr e 



-/W/(7.(*)) 



(3.15) 



using the rescaled Hamiltonian i?'({cr-}, -f[ , 7 2 , N') and 
the H'({ai}, 71, 72, N). Since the Hamiltonian is of the 
same form as that before the scale transformation, so 
does the free energy. Comparing the coefficient function 
of the spin-spin coupling term of the effective Hamil- 
tonian, we get the transformation function Ul- In the 
vicinity of the non-trivial fixed point 7* , 
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(3.16) 



we perform Taylor expansion around 7* = (K*, 72), and 
make a truncation to the first order (the simplest case). 
The renormalization group transformation is identical 
with coordination transformation, 



71 - 7i 

72 - 72* 



(pk) (pk) 

v a 7 i ' v o 7 2 ' 



7i 
72 



7i 

72 



(3.17) 



We denote 87' — (7^ — 7*,7 2 — 7 2 ) T and £7 = (71 — 
7i, 72 — 72 ) T , the element of the group is 



/ (M) (M) 

v 7 1 ' v 9 7 2 ' 



(3.18) 



This is the first order approximation of the game opera- 
tor, the exact game operator S-y' = U8y may be obtained 
by including the higher order approximations. 

If we take this Ising model as a war game, the full gen- 
erals of the two army are 71 and 72 . The spins confined in 
the lattice sites are soldiers. The two generals take better 
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and better strategies to play through scale transforma- 
tions. The player delivered his message to his opponent 
through scale transformation matrix. Then they adjust 
their physical parameters 71 and 72 in the next round 
of game. This game process is represented by a series 
of game operator, /, Ul, {Ul) 2 , {Ul) 3 , ■ • • • The game op- 
erator actually defines a flow from high energy to low 
energy by the change of scale. 

From physicist's point of view, high energy means 
hight momentum. While the momentum is characterized 
by the fourier transformation of lattice spacing on trian- 
gular lattice. If we divide the lattice space smaller and 
smaller, its dominant momentum representation grows 
higher and higher, and finally leads to divergency in the 
continuum limit. In the war game, the full general may 
roughly divide his army into tens of corps, and he is only 
in charge of the tens of major generals, this is the low en- 
ergy case for a physicist. The full general may continue 
to dived his army into hundred of regiments, and further 
into thousands of companies. If he is powerful enough 
like god, he can directly take charge of the millions of 
soldiers, this is the high energy part of a war game. 

Effective low-energy theories can always be reached by 
integrating the high-energy degrees of freedom. This is 
not merely a conception in physics, it occurs every day 
everywhere in economics . If we went to free market, 
the bargaining procedure between sellers and buyers is 
actually a dynamical demonstration of renormalization 
group transformation. In the beginning, the seller would 
present a very high price to maximize his profit, the 
buyer would try to lower it down in order to decrease 
his damage. This is what happens: Buyer:" How much?" , 
Seller:" 1000 dollars." , 

Buyer:" Too expensive, if you sell it at 500 dollars, I would 
buy it.", 

Seller:" No way, how about 500+±500 dollars", 
Buyer:" 500+^500-^500!", 
Seller:" 500+|500-i500+|500!", 
Buyer:" 500+±500-±500+±500-^500!", 

Buyer:" 500 (1+lim^oo E»=i M)" +1 2-")!" , 
Seller:" Done!". 

The fixed point of this bargain is done when N — > 00. 
Usually this bargain does not go so far, people always 
cut it off after two or three round of bargains. Here we 
made the assumption that both the buyer and seller are 
rational. In reality, there are various cases. If both the 
buyer and seller are adamant, they may trapped in a 
p-period circular bargain, i.e., B: "1000 dollars", S:"500 
dollars", B: " 1000 dollars" , S: "500 dollars" ,• • • , B: "1000 
dollars", S:"500 dollars", • ■ ■ . In a more complex case, 
the buyer bargains following a complex mapping rule, 
his pth offer is 7° = (7^ + 7 b p _ 1 + j b p _ 2 + ■■■), where 

7p_ 1 is the seller's (p — l)th offer, the game may end up 
with chaos, bifurcation, all kinds of nonlinear phenomena 
come in. In fact, these nonlinear interacting phenomena 
indeed occur during a phase transition. 



The renormalization group transformation is just the 
bargain game rule. The bargainers in physical system are 
different interactions. Practical physical system usually 
has good negotiators. The bargain series converged into 
a fixed point. The unit money on which the buyer and 
seller are bargaining becomes smaller and smaller dur- 
ing the dynamic process of renormalization group trans- 
formation. In physics, this means more and more high 
energy effects are integrate into the effective low energy 
theory. 

Attractive fixed points describe stable phases within 
the renormalization group formalism. The critical point 
of phase transition corresponds to saddle point at which 
the physical parameters in the game reach equilibrium, 
namely the Nash equilibrium. At the Nash equilibrium 
point, if any one of the player take a wrong strategy, he 
would fail and flows to the attractive fixed point, then 
he is confined in a stable phase. The Nash equilibrium 
point is where all different phases coexist. 

There is a theorem in game theory concerning on 
the existence of the saddle fixed point. Applying 
the Brouwer's fixed-point theorem, Aubin obtains a 
corollary Q: Suppose that the behaviors of 71 and 72 are 
described by one-to-one continuous decision rules and 
that the strategy set {71} and {72} are convex compact 
subsets of finite-dimensional vector space, then there 
is at least one consistent pair, i.e., there exist at least 
one Nash equilibrium. This theorem may help us to see 
whether phase transition exist or not for a given game. 

The application of Renormalization group to game 
theory an example: Cournot duopoly 

Usually it is very hard to find the Nash equilib- 
rium of multi-player games. Renormalization group 
theory provide us new tools to find the Nash equilibrium 
solution for a multi-player game. We first transfer the 
multi-player to an effective quantum or classical many 
body system, and find its Hamiltonian, we can apply 
various well-developed numerical renormalization group 
calculation method to find the critical point of phase 
transition, then we find the Nash equilibrium solution. 

We consider Cournot duopoly game (see Appendix H). 
The simplest case is there are only two players Alice and 
Bob, they are manufacturers of the same kind of product. 
In the beginning, Alice produces 7$, and Bob produces 
7q. The two payers competes with each other, each of 
them changes their productions according to their oppo- 
nent's production. One chooses proper strategy to min- 
imize his own cost function. Alice's canonical decision 
rule is 7^ = /a(7o)' 7i = /b(7o)j we express this game 
into matrix form 

(3)-tt*)(D- 

The second round game follows 7^ = JU(7i) = 
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/a(/b(7o°)). ll = hilt) = fB(f A (l b )), i- 
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7o 6 



(3.20) 



This equation has a more clear expression 72 = U 2 Jo 
after we introduced the game operator U and strategy 
vector 7, 



[7 



f A 

h 



7« 



7„ 

7 6 



(3.21) 



Suppose the players paly alternatively, Alice is in the 
even periods and Bob in the odd periods, when Bob 
produces 72„_i in the period 2n — 1, Alice produces 
72« = /a(72«-i) m the period 2n, Alice then changes 
her production rate and produces 72„+i = /s(72n)i an d 
so on. The sequences of 7f n and "i\n-i satisfies the re- 
cursion relation, 



(3.22) 



This sequences converge to |. ^ is the fixed point of 
U n+1 ,(n^> 00), 



lim 7„ = lim ^" +1 7o =7*. 



(3.23) 



Now we see this game operator is the same as bargain 
game in the sense of the renormalization group transfor- 
mation. 



3.2. Solutions of Renormalization group 
transformation equation in game theory 

The game in reality always experience a dynamic pro- 
cess. For example, the bargain on price between the seller 
and buyer still exist in primary market. The advent of su- 
permarket drive this kind of observable negotiation pro- 
cess to the backstage market investigation. The operator 
of the supermarket adjust its next day's merchandize dis- 
tribution according to last day's sell. A good customer 
also compares today's price with previous price to buy 
what he needs at a lower price. This is a mutual inter- 
acting process which may be expressed by a nonlinear 
game operator, 



7„ 



f A 

h 



7," 



(3.24) 



where 7 is the seller, 7 & is the buyer. This game has an 
equivalent representation by a pair of difference equation, 



In = fAhLl), 



(3.25) 



In fact, the buyer not only consider the seller's price, 
but also consider money amount he could pay. On the 
other hand, the seller must reflet the amount and quality 



of the product or service he could offer. So a complete 
difference equations of this game is 

In = /A(7^-l, In-ll In = fs(^-l, 7„-l)- (3-26) 

In fact, the players of a game mainly concerns about 
profit difference between last round and next round, the 
most efficient way is to take the profit of last round as a 
datum mark. So we can always find a term like y%_i on 
the right hand of Eq. (|3.26p . Thus the game difference 
equation is transformed into differential equations, 



/a 

In 



7 



it _ u 

In In 



lim 

ft-vo h 



lim ■ 

h~*0 



l b n-l 



^fl(7«-i»7»-i), (3-27) 



h is the step size of the tuning parameter. This equation 
may be viewed as renormalization group transformation 
equation. The two functions Fa/b form the game op- 
erator which map one state in the strategy space into 
another. Its trajectory depicts the renormalization flow. 
The strategy space of two player game is two dimensional 
parameter space. The corresponding game operator is a 
general 2x2 matrix. This two dimensional matrix is an 
element of renormalization group. For some special case, 
if the game operator is an element of spin(3), we can 
expand it in terms of traceless Pauli matrix. 

The game operator of a iV-player game is a traceless 
N x N matrix, 



U = 



( hi 

hi 



f } 2 
/22 



J2N 



(3.28) 



V /jVl fN2 ■ ■ ■ fNN I 



Following the same procedure as two player game, we can 
get an equivalent N dimensional difference equations. We 
summarize the output field and input field into one state 
vector 7 = (O u t,fi n ), a general system is governed by 
the differential equations, 



9 s 7 l = L(7, <9 7 )7i 



1,2,. 



(3.29) 



where s is the tuning parameter, it could be any input 
field. In fact, the output and input are not absolutely dis- 
tinguished, that depends on what we have known, what 
we still do not know. We usually take those we can ma- 
nipulate as input, and those we will detect as output. 
When physicist study phase transition of a system, they 
usually tune some external parameter, such as temper- 
ature, magnetic field, et al. They investigate how the 
system response according to different physical parame- 
ters so that we may control it for practical purpose. For 
example, superconductor physicist study how the con- 
ductivity behaves when people raise the temperature as 
high as room temperature, they do not care much about 
the time. If the parameter s is taken as time, it is just 
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the conventional dynamic system, a physicist's a ppr oach 
to its solution is the algebra dynamic arithmetic [10(. 

The game operator L(-f, <9 7 ) is the infinitesimal gener- 
ator of the translational transformation with respect to 
parameter s. For a given initial strategy 7 , the evolution 
of the system along s follows 

(s) = e H^) s = J2 £!l«( 7; a 7 )/. (3.30) 
i n! 

This equation enveloped the whole process of renormal- 
ization group transformations. Take the the bargain 
game as an example, n = is the first round of bargain, 
the buyer and seller put their initial cards on the table. 
n = 1 is the second round of bargain, after they have set- 
tled the main business down, they began the bargaining 
on minor business now, the amplitude of this round is 1. 
The amplitude of the third round of bargain is 1/2!. As 
the renormalization group transformation goes on, the 
amplitude of the nth round of bargain decays following 
1/nl. When they made an agrement on all the problems 
from the dominant ones to the ignorable ones, peace ar- 
rived. 

The ending point of this game is fixed point which is 
given by the exact evolution state vector 7(s). The un- 
stable fixed point corresponds to the Nash equilibrium 
point, any one of the players takes one more aggressive 
step to increase his own profit would result in the break- 
down of the deal. The Nash equilibrium point is where 
the phase transition occurs. The stable fixed point indi- 
cates stable phase. 

3.3. Symmetry losing as a classification of phase 
transition 

We gave a very general definition of phase transition 
in section ()2.2|) . phase transition is a game between old 
phase and new phase. Whenever the winner becomes 
loser or vices versa, phase transition occurs. Transition 
is always accompanied by the transfer of award from loser 
to winner. The award is quantized, so we observe sudden 
change of output across the phase boundary. 

Each time we make a renormalization group transfor- 
mation, the players accomplished one round of game. If 
someone lose, someone must win. The temporary phase 
boundary invades from the winner into the loser. The 
amplitude of boundary's change in the nth round of game 
is proportional to 1/n!. If n = 1, that is the 1th order 
phase transition, the amplitude of sudden change is most 
significant, the phase boundary is roughly fixed except 
some unsettled regions. Then we need the second round 
of combat to negotiate the main part of the unsettled re- 
gion. As this kind of negotiation goes on to higher order, 
the unsettled boundary becomes smaller and smaller, its 
amplitude decays following 1/n!. When all the unsettled 
boundary are fixed, we reached a fixed point. If this fixed 
point is stable, we are in a stable phase. If this is an un- 
stable fixed point, we are at a critical state, to be or not 



to be, it lies in the hand of this point. Any minor devia- 
tion from this point would decide who is winner, who is 
loser. 

In the region of stable phase, the state vector 7 of the 
game is C p continuous function, i.e., the derivative of the 
state vector -r-£ is continuous up to the pth order. This 
differentiability of the pth order derivative beaks at some 
singular points, at which we would observe a pth order 
phase transition. 

We usually take the pth order non-differentiability of 
output vector to measure phase transition. These output 
vectors O u t(u) are the subset of the state vector, usually 
they are the cost function of the game. The control vec- 
tors \u) are game players. For a quantum many body 
system, the output vectors O u t(u) could be any statis- 
tical observables or any external response, such as free 
energy 0\ — F, ground state energy O2 = E g , thermal 
potential O3 = Q, susceptibility (9 4 = x, specific heat 
O5 = Ch, correlation length Oq = £, compressibility 
O7 = kt, ■ ■ ■ ■ The players of the game is control vector 
its component includes all input variables, for ex- 
ample, temperature u\ — T, pressure u-i — P, effective 
external magnetic U3 = 72 = fiB/ksT, spin-spin inter- 
action Ui = 71 = J/ksT, electric field U5 = E, chemical 
potential U5 = n,- • • . 

The Nash equilibrium of this A^-player game is to find 
the eigenvector of the operator U n+1 in strategy space. 
This is tantamount to solve the equation 

(U n+1 -I)\ip) =0. (3.31) 

This equation indicates an infinitesimal transformation 
around the identity. A finite transformation is con- 
structed by the repeated application of this infinitesimal 
transformation. For a realistic game, it is always impos- 
sible to get an exact matrix of linin^oo U n+1 up to the 
infinite order. The output function encounter divergence 
on the input base manifold. These singular points are 
where the transition from loser to winner occurs. Out of 
these singular regions, the output function has very good 
behavior. 

Renormalization group is a semigroup, because its el- 
ements have no inverse. The ordinary group is a subset 
of the Renormalization group. So they only provide us 
some basic understanding to phase transition in certain 
special cases. 

In critical phenomena, the correlation length between 
particles goes to infinity at the phase transition point 0], 
there is an obvious conformal transformation symmetry. 
This is a kind of symmetry when we take the physical 
particles as players of a game. But a practical physical 
system always contain billions of particles, this is not a 
good choice. 

A more practical way to study quantum many body 
system as a game is taking the different interactions as 
players. The interaction parameters form the strategy 
vector of the game. As we know, the strategy vector 
at fixed point is invariant under the operation of game 
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operator. The holonomy group at the Nash equilibrium 
point is a special subset of the game operator space. We 
may check the discontinuity of the output field under 
holonomy group to check the phase transition point. The 
strategy space is expanded by the strategies of all play- 
ers, the holonomy group transformation is actually the 
transformation on the strategy manifold. 

Lie group is more familiar to most physicists, it is also 
a subset of renormalization group. When the game has 
a continuum of players, the strategy space construct a 
manifold. The payoff functions is a cross section of the 
fibre bundle established on this base manifold. We can 
define the order of phase transition from the loss of Lie 
group symmetry 

We take the two dimensional Ising model as an exam- 
ple to demonstrate the explicit procedure. The output 
vector could be chosen to have only one component — the 
free energy, and the input vectors are 71 and 72 , they are 
correspondingly the coupling interaction and magnetic 
field. 71 and 72 choose the best strategy at each step to 
increase his own welfare and decrease his loss. They first 
choose the initial strategy pair of arbitrary value. Then 
they behaves following the decision matrix in the next 
round of game. Repeating n rounds of this game, they 
reach the fixed point. At the Nash equilibrium point, 
both the two parameters has no further steps to increase 
his welfare any more. 

There are only two independent parameters under the 
constrain of the equation of state, any finite transforma- 
tion around the fixed point could be reached by repeating 
infinitesimal of 5*0(2) whose generator is 



f d d 

L = 72^ 71-5—- 

071 C>72 



(3.32) 

The group element of SO(2) is expressed by the expo- 
nential map, Lie : I^>Exp(I), 



U(6) 



„9L 



= I + L6 + -L 2 9 2 



(3.33) 



We expand the group element up to the pth order, i.e, 
U(0) = Eo^(^ 6 P ' the definition for the pth order of 
phase transition reads, 

p 



(I 



j2l(L9rd^Y,-(L6rd 



(3.34) 



As p — ► 00, it reaches the exact 5*0(2) group element 
U{6). When we choose the output vector as free energy 
F, and the input vector 71 as temperature T and 72 as 
pressure P, this definition of phase transition has unified 
all orders of Ehrenfest's definition into one equation. For 



example, for n — 0, — I, it yields F A 
n = 1, U = I + L6, then (I + L9)F A = (I + 
derived 



= F B . For 

L9)F B , 



we 



=> F A = F B , 



Q F B 

dT 



dF A dF B 

= 0, 



dT 



dP 



dF A 
~~dP 



= 0. (3.35) 



Therefore the essence of Ehrenfest's definition for dif- 
ferent order of phase transition actually depend on how 
many order of the symmetry of free energy is pre- 
served during the phase transition. We denote — 
YZ ji(L6) p , Eq, ((334)) may be reduced to 

U^O A t = U^O B , U^6 A ^U^6 B t , (3.36) 

So when p — > 00, it reaches the exact 50(2) transforma- 
tion. That means the output vector is differentialable to 
infinite order. Since there is no discontinuity, we are not 
able to observe it from external responses. 

Eq. (|3.36p has provided us qualitative understanding 
to how the symmetry loss induced a phase transition. It 
hold for the more general case that the output vector is a 
cross section of fibre bundle on a manifold expanded by 
many parameters. For example, if O u t is a vector field of 
(71, 72,73), the simplest choice is to introduce the gener- 
ators of 50(3) which are the three angular momentum 
operator, 

d d 
M = -73^ h 72-^—, 

072 073 

d d 
^2 = -71-5 1" 73^—, 

073 C71 

L 3 = -li-K— + 7i J—, ( 3 - 37 ) 
071 072 

The order of phase transition is characterized by the 
group element of 50(3), U = e w "' L . We expand U to 
the pth order, and investigate the equation 



5*>[U(0)(Oi 



Out)} 



56p 



9=0 



SP[U(0)SOut] 
S0P 



8=0- 



(3.38) 



If the field tj> ^ for the pth order transformation, but 
vanishes for all the transformation under the pth order, 
we define the phase transition as the pth order. 

The game operator L in renormalization group ele- 
ments U = e l6n ' L could take arbitrary sophisticate for- 
mulas. L can be expressed as polynomial operator ex- 
panded by 7 and 9 7 , 



0'' 



(3.39) 



The specific form of L relies on specific systems. No mat- 
ter how complex it is, we can always check the expansion 
of the exact solution 7(5) = e i( 7> <9 -') s 7 to fond out the 
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singular points which separate the whole space into dis- 
crete regions. 

It must be pointed out here, the symmetry loss here 
is different from the conventional spontaneous symmetry 
breaking in physics. Take the familiar ip A model for ex- 
ample, there is a spontaneous symmetry breaking from 
SU(2) to U(l) for vacuum state, it is the lagrangian of 
tP 4 model C = ±(<9 M V) 2 - |mV " \^ that has SU(2) 
symmetry. But when the ip 4 model is studied using the 
game theory of renormalization group transformation de- 
veloped in this paper, we do no care about the i/>-field at 
all. We just take the mass m and coupling constants A 
as two players, and take the physical observables calcu- 
lated from the partition function as output vector. It is 
on the manifold expanded by the mass and coupling con- 
stant we introduce the SO (2) transformation around the 
critical point. 



4. THE TOPOLOGICAL THEORY OF 
UNIVERSAL PHASE TRANSITION 

In physics, the renormalization group transformation 
is going on under the constrain that the Hamiltonian of 
the system must has the same form after transforma- 
tion as it is before transformation at the critical point, 
= ^(7i)72)- So the partition function and 
free energy function also maintain their form during the 
transformation. Following the spirit of special relativity, 
Einstein's principle of general covariance states that all 
coordinate systems are equivalent for the formulation of 
the general laws of nature. Mathematically, this suggests 
that the laws of physics should be tensor equations. In 
this sense, the laws that governs the motion of everything 
in universe should not depend on coordination, no matter 
it is in physical world or social world. 

Phase transition is perhaps the most common phenom- 
ena in nature as well as in human society. The basic law 
of phase transition does not depend on its base manifold, 
on which we established the equation of the states for a 
certain system, i.e., physical system, chemistry system, 
biology system, or social system. What we face is a black 
box, the only source of information about the inside of 
the black box is the output vector which responses when 
we alternate the input. 

When we confine our issue in physical system, the out- 
put vectors O u t(u) are macroscopic observables which 
may be decided in the frame work of physical science or 
directly measured by conducting experiment, such as free 
energy 0\ — F, ground state energy 2 = Eg, thermal 
potential O3 = fi, susceptibility O4 = x, specific heat 
O5 = Ch, correlation length 6 = compressibility 
Oj = kt, ••• ■ The input vector are also macroscopic ob- 
servables, such as temperature R\ — T, pressure R2 = P. 
particle density R3 — N, volume R4 = V, chemical po- 
tential i? 5 = /v • • , conductivity Rj = a, susceptibility 
Rk = Xi an d so on - The input vector and output vector 
are relative, they are no fixed once for ever, it depends on 



our subject. For instance, if we study how the volume of 
a gas changes when we change the pressure, the volume 
is the output vector, and the pressure is the input. If we 
intend to study the inverse relation between volume and 
pressure, the input and the output exchange their roles. 

In this section, we choose the base manifold of output 
vectors. Our interest focus on the intrinsic geometric and 
topological quantities, for they do not rely on the local 
coordination system. 



4.1. Topological current theory of phase transition 

The input vectors are players of a game, the output 
vectors are the cost function or payoff function. The in- 
puts take different values to maximize their profits, the 
order of the strategy they played during the game is of 
crucial importance. About 2300 years ago in Chinese his- 
tory, general Sun play horse racing with the King, both of 
the two players have three horses, a weak one, a regular 
one and a strong one. A match consists of three rounds, 
each of the three horses must take part in at least one 
round. The King's horse in each class is more power- 
ful than that of Sun's. But Sun finally won, he used 
his regular horse race against the King's strong horse, 
he certainly lost this round. But his regular beat the 
King's weak, and his strong beat the King's regular. If 
Sun change the order of any two of his horses, he would 
lose the game. Different order of strategy lead to totally 
different results. 

A physical process always depends on many external 
parameters which entangled with one another. These ex- 
ternal input are game players, their different values rep- 
resent different strategy. If we try to liquify a realistic 
gas confined in a rubber container by applying high pres- 
sure and cooling. There are two ways, one way is first 
to press it, and then to cool it down, the other way is 
first cooling it and then pressing it. These two ways do 
not always lead to the exact same structure except some 
special cases. More over, the speed of the cooling has in- 
evitable effect to the final result. A medium carbon steel 
is transformed to austenite at about 1,550 degrees. If it is 
allowed to cool slowly back to room temperature, it has 
the ferritic structure. If it is rapidly cooled, the austen- 
ite quenched to another shape which has high-strength 
structure. 

So a general output vector field is established on curved 
manifold, for physicist, the output vectors are physical 
observables on curved parameter space. The strategies 
are not commutable in the curved output space. The 
physical parameters are commutable only if the base 
manifold is homeomorphic to flat Euclidean space. Its 
physical representation is a description of adiabatic phys- 
ical transition process which is reversible in classical ther- 
modynamics. 

To study the intrinsic geometric properties of phase 
transition, we choose the most fundamental manifold of 
output field O ut (7)- The flow of vector field of renormal- 
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ization group transformation point out where the Nash 
equilibrium point is. If the vector field is continuous ev- 
erywhere on the whole base manifold, one observe noth- 
ing. If only there appears a discontinuity, then we make 
sure that there is something happening. The fundamen- 
tal vector field to detect the discontinuity is 



50p 



(4.1) 



We first take two component of this vector field to study 
the pth order phase transition. They are the pth order 
tangent vector field on O u t("f), 
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O ut (7), 



(4.2) 



where p is arbitrary number, jj are input vectors. Here 
we first take the simplest case that has only two input 
field as an example to present the basic relation between 
topology and phase transition. 

In the neighborhood of the critical point, the manifold 
is approximately holomorphic to flat Euclidean space, the 
translation operator <9 7 is good enough to describe the 
transportation of output field. But topology concerns 
about the global geometry of the manifold, we have to 
walk out of the vicinity of the saddle point, then we need 
to introduce the covariant derivative, 



n — r) 



D~f 2 — ^72 *^72 ' 



(4.3) 



A lt and A 72 is the gauge potential which connects the 
vector field between different domains. The covariant 
derivative of the tangent vector field on O u t (7) is 
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(4.4) 



The commutator of the two covariant derivative produce 
the Gaussian curvature on the two dimensional Rieman- 
nian manifold, [Di,Dj] = — iQy with where Qij = dij — 
diAj — djAi. The Gaussian curvature may decompose 
into the product of two principal curvatures, «i(7 ) and 
^2(7°). When m > and K2 > 0, it is a elliptic surface, 
when «i > and K2 = 0, it is parabolic, for n\ > and 
K2 < 0, it is hyperbolic. If the Gaussian curvature is 
zero, it means the manifold is flat. 

As it is well known, Euler characteristic number on the 
compact two-dimensional surface is defined by Gaussian 
curvature G, i.e., — f M G^/gcPx. In the differ- 

ential geometry, the Euler characteristic is just a special 
case of the Gauss-Bonnet-Chern theorem which defines 
a topological invariant on the 2n dimensional surface. In 
terms of the Riemannian curvature tensor, the Gaussian 
curvature fl is written as 



Cl = -— -pRfiuXa- 



(4.5) 



As all know, the Riemannian curvature tensor correspond 
to the gauge field tensor in gauge field theory [llj. When 



we set up the vielbein e° on the surface, the Riemannian 
curvature tensor can be expressed in terms of the gauge 
field tensor as R^xa = -e^e*F°J. We introduce the 
50(2) spin connection for the tangent vector Eq. 
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(4.6) 



Then we can define an Gauss mapping n with n a — 
n a n a = 1, where \\(j)\\ = y/p^ (a = 1,2). 
Considering the symmetry of the unit vector field n, we 
introduce the SO (2) spin connection the covariant 
derivative of n is defined as D^n a = d^n a — uj^ b n b . Notic- 
ing SO(2) is homeomorphic to U(l). There is a one to 
one correspondence between SO (2) spin connection and 
U(l) connection. Both of them have only one indepen- 
dent component. As shown by Duan et al[l2j|. one con- 
sider parallel transportation from Di<f> = di4> — iAi<p = 
and its conjugate equation = di<ft +iAi<f>^ = 0, it is 
easy to obtain the U(l) gauge potential Ai = e a bn a din b . 

It was proved that the Gaussian curvature G can be 
expressed as 



CI 



1 e" w 



F tlu = (d fl A u -d u A IJ ,), (4.7) 



here F^ u is U(l) gauge field tensor. Substitute the U(l) 
gauge potential Ai = € ao n a din b into the Gaussian cur- 
vature, one may express the Gaussian curvature into a 
topological current 



0= £ 

i,j,a,b—l 



e tJ e ab 



dn a dn b 



(4.8) 



This topological current appeared in a lot of condensed 
matter systems. In differential geometry, the integral of 
the gauge field 2-form is the first Chern number G\ — 
J M fl, it is the Euler characteristic number on a compact 
Riemannian manifold. 

From the unit vector field n, one sees that the zero 
points of </> are the singular points of the unit vec- 
tor field n which describes a 1-sphere in the <fi vec- 
tor space. In light of Duan's <fi— mapping topological 
current theory fH), using = ^ + <t> a di-^ and 

the Green function relation in </>— space : 9 a 9 a ln||0|| = 
27r<5 2 (0), (d a — j^), one can prove that 

Vl = 5 2 {$)D{t) = 5^$){ (t) \^}, (4.9) 
7 

where D{<f>/q) = \ Y.l 3 .a.b=i^ k ^bd^ a d k ^ b is the Ja- 
cobian vector. In the extra-two dimensional space, this 
Jacobian vector is just the Poisson bracket of (j) 1 and <p 2 , 
W\0 2 } = Ei($£^ - As we defined above, 

the vector field <f> is the tangent vector field on the out- 
put manifold. This tangent vector field may be viewed 
as a projection of a source vector field, which satisfy 
^ cos 9 = (f>, 9 is angle between the vector and the tan- 
gent plane at point p, Obviously when <f> — 0, the source 
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vector field points vertically up. If we draw the configu- 
ration of the vector field around point p, one would see 
a Skyrmion configuration. There is a sharp peak around 
the point at which cf> — 0. These sharp peaks bears a 
topological origin. 

Eq. (|4.9p provides us an important conclusion imme- 
diately: n = 0, iff $ ^ 0; Q i= 0, iff = 0. In other 
words, the solutions of the equations 

0! = ^ 1 - 1 - Q ^- 1 + a O„ t (7) = 0, (a ± b) 

02 = d^-Od^+'Outtf) = 0, (4.10) 

decides whether the phase transition exist or not. If 
Eq. (|4.10p has solutions, it means that there are some 
points on the domain wall at which the tangent vector 
field is continuous. According to our general definition 
about phase transition, a phase transition is a revolu- 
tion, it happens when the system can not survival with- 
out changing itself to fit the new environment. If there is 
any strategy that the system could take to survive, it will 
not take any risk to face revolution. Therefore as long 
as there exist solutions for = 0, we will not observe 
any sudden change, this state is marked by a non-zero 
Euler number. On the contrary, if Eq. (|4. 10[) has no 
solution over all strategy space, this means the tangent 
vector field across the domain wall encounter a barrier, 
it has to jump over the barrier with finite hight to access 
another domain. In the game theory, this indicates the 
system can not find any strategy to help itself move from 
one domain to another domain, a reform or revolution is 
required. This indicates a phase transition, this state of 
system is marked by a zero Euler number. 

Now we see Eq. (|4.9|) actually describes topological 
configuration of vector flow around the surviving strat- 
egy. Each surviving strategy present a peak on the tan- 
gent vector plane. More peaks means longer life time for 
the old phase, if peaks become less and less, that means 
the old phase is dying. The system becomes more and 
more unstable, and began to collapse. When there is no 
peak left, the system has to start a revolution up, the 
system is totally unstable now, chaos effect come into 
action. 

The number of solutions of equation = counts 
the number of surviving strategy for the old phase. 
The implicit function theory shows, under the regular 
condition [ 13| D(0/7) ^ 0, we can solve the equations 
= and derive n isolated solutions, which is denoted 
as Zk — (7fe,7fc), (k — 1, 2, ... n) At the critical point Zk, 
the Jacobian D(^) can be expressed by Hessian matrix 

M <5o„ t (7) of <50„ t (7*), i.e., D{^)\ Zk = detM SOut (i){zk)- 
According to the 6— function theory [lj], one can expand 
5(0) at these solutions, <5 2 (0) = ^Li Pk tSkP ■ Then 
the topological current becomes 



i 

k=l 



detM SOut(i) {z k ) 2 
Pk \detM SOut(i) (z k ) r 7 



here pk is the Hopf index and the Brouwer degree 
i] k =signD(0/g) 2fc =±1, i.e., signdetM &F (z k ) = ±1. 
.Applying Morse theory, we can obtain the topological 
charge of the transition points from Eq. (|4.8[) . 

Following Duan's 0— mapping method, it is easy to 
prove that PkVk — W k is the winding number around 
the kth. critical point. The winding number measure how 
many times the vector flow surround the isolated surviv- 
ing strategy. When the output manifold extend a com- 
pact oriented Ricmannian manifold in strategy space, the 
total winding number 



Ch = / VLd 2 q = J2 Wk 



(4.12) 



k=l 



(4.11) 



is just the Euler number. Euler number is a topolog- 
ical number, it strongly relies on the topology of the 
base manifold. Especially when the output manifold is a 
compact orientable Riemannian manifold, such as sphere, 
torus, or a disc with boundary, and so on, the total topo- 
logical charge of the transition points is the Euler num- 
ber. The Euler number of a 2-sphere is Ch — 2. Accord- 
ing to Eq. (|4.12p . we see if there are two peaks of the 
vector field distributed on output manifold, each point is 
assigned with a winding number W = 1 to make sure the 
Euler characteristic number of the 2-sphere. If there is 
only one peak, the winding number must be 2. 

As we know, each peak represents a surviving strategy, 
more strategies provide more surviving opportunities for 
the old phase. But there is a topological constrain from 
the base manifold which requires that the sum of winding 
number around each strategy must be equal to the Euler 
number. If we require the winding number must be posi- 
tive, we see if there are four strategies, each of them must 
carries half winding number Wi = 1/2, (i = 1, 2, 3,4). 
The optimal distribution of the four strategies should 
make them separated as far as possible. For if they are 
at a crowd, it would be dangerous, their enemy does not 
need to spend much energy to block them all, then the 
old phase dies. The four critical points of the same sign 
repel each other to reach the minimal of the system's 
total energy, so they will be separated as far as possi- 
ble. When the equilibrium is reached, the most likely 
distribution is the four critical points are situated at the 
vertices of a tetrahedron. As for the two strategies case, 
one sits at the North Pole, the other sits at the South 
Pole. If there is only one strategy with W = 2, it is 
unstable and is apt to split into two or four. If there 
is a strategy with W = +3, the strategy with negative 
winding number would appear, but these state are very 
unstable. 

The output manifold may jump from a sphere to a 
torus, and to a torus with many holes, the Euler num- 
ber would jumps from Ch = 2 to Ch = 0, then to 
Ch = 2(1 — h) with h as the number of holes of the 
torus. The topological change of base manifold would 
either kill the old phase or save its life, so topology plays 
a very important role in phase transition. 
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(a) 



(b) 



FIG. 2: (a) The distribution of four surviving strategy with 
topological charge +1/2 on a 2-sphere. (b) Two surviving 
strategies with topological charge W — +1. 



What we presented in the discussions above is the sim- 
plest case, there is only one output field with two input 
vectors. For the most general case, the output vector has 
to components with to input vectors. There exist a to 
dimensional tangent vector field for every component of 
the output vector, 



where X)j a » = 0- We can define a gauss map, i.e., an 
unit vector field 77 with n a = a /||</>||, n a 77 a = 1, where 
1 10| I = \/4) a 6 a (a — 1, 2, 777). Following the topological 
field theory [l2j, we can find a topological current, 



CI = 



E 



ij...k ^ 

e e a b.. 



,c=l 



977 a dn b dn c 



(4.13) 



where (a,b,...,c = 1, 2, 777), (i,j,...,k = 1, 2, m), 
£ab...c is antisymmetric tensor. On even dimensional 
manifold, this topological current is exactly equivalent 
to the Riemanian curvature tensor which directly leads 
to the Gaussian curvature in two dimensions. Applying 
Laplacian Green function relation, it can be proved that 
f2 = 5((f>)D(-), where the Jacobian D(&) is defined as 



D{~) 

7 



E 



z,j,...,fc;a,fc,...,c— 1 



In to = 2 71 dimensional manifold, it was proved that the 
topological charge of this current is the Chern number 
Ch = J Vld 2 q = Y^,k=i Wk, which is the sum of the wind- 
ing number around the surviving strategies for multi- 
player game. 



4.2. The universal equation of coexistence curve in 
phase diagram 

A phase transition is a war, is a game, is a revolution. 
No matter where it takes place, it becomes landmark in 



history of a system. The renormalization group trans- 
formation theory told us a phase transition point is the 
Nash equilibrium solution of a game. As shown in the 
bargain game, the profit of seller is the loss of the buyer 
and vice versa. When the seller takes proper strategies 
to maximize his profit, the buyer is trying to minimize it 
by taking strategies from a different space. So the seller 
is approaching to the maximal point of the payoff func- 
tion, in the meantime the buyer is looking for its minimal 
point. A Nash equilibrium appears at their intersection. 
The Nash equilibrium point is a saddle point, the output 
field reaches its maximal point in 71 direction, but get a 
minimal value in 72 direction. The derivative of the out- 
put field corresponds to 71 and 72 must be of opposite 
sign. This leads to the coexistence equation for different 
phases. 

In previous sections, when solving the equation of tan- 
gent vector field <j> = to find the surviving strategies, 
we applied a regular condition D(6[j) ^ 0, which comes 
from the implicit function theorem[13l|. When the regular 
condition is violated, i.e., D(<fi/q) = 0, a definite solution 
of equation <f> = is not available. Then the branch pro- 
cess of the solutions function occurs. A mathematical 
demonstration of the branch process could be found in 
Ref. [13] • This bifurcation may be understand from game 
theory. Under the regular condition D(cf)/j) 7^ 0, if the 
solutions of <fi = exist, that means the old phase still 
has strategies to survive, if there is no solutions, the old 
phase can not find any strategy to make a living, it has 
to die. When the regular condition fails, D(<fi/j) = 0, 
the old phase and new phase is at an equilibrium war, 
if the old phase win, it find ways to survive, if it loses, 
no surviving strategy exist, the old phase dies. There- 
fore, it is at the very battlefield of D(cj)/^) = 0, the two 
phases have equivalent power, nobody wins, nobody lose, 
but they are fighting against each other. Several roads 
branched out of this battlefield, to be or not to be, the 
old phase has to make a choice when passing this critical 
region. 

In two dimensional input space, the coexistence 
curve equation D((p/j) — is just the familiar 
Poisson bracket for the tangent vector field 0* = 



-1— ii . 



-l-ia 



{0\0 2 } = O. 



(4.14) 



It is an unification of the special coexistence equations of 
different order phase transition. As all know, in quantum 
mechanics, if two operator A and B commutate with each 
other, i.e., [A, B] = 0, they share the same eigenfunction. 
Eq. (|4.14p means that the two classical field (j) 1 and (j) 2 
are commutable at the phase transition point. 

When the game has three players, the output field is 
a function of three parameters, each of them holds the 
life of an old phase. The three phases intersects with 
one another at the coexistence points which sit at the 
solutions of 



{0 l ,^,0 fc } = O, 



(4.15) 
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where {(/>*, J , <j> k } is the generalized Poisson bracket. Its 
quantum correspondence id the Jacobi identity 



[A, [B,C\] + [C, [A,B]] + [B,[C,A]]=Q. 



(4.16) 



For a n- player game, we need to introduce a n- 
dimensional renormalization group transformation on the 
output manifold. The transformation operator expand 
the tangent vector space around the identity on the man- 
ifold. We denote a vector operator as L, a group element 
is given by U — e i6n ' L . The basic tangent vector field for 
phase transition is 



y 7l , V72 1 ' 



SP[U(9)S(0\6 \0}} 
60P 



9=0- 



(4.17) 



The most general definition of generalized Poisson 
bracket [l5| for n component vector field is 



{4> 1 ,<t> 2 ,...,<j> n } = 



d^\4>\-,4> n ) 

9(71,72, -,7n) ' 



(4.18) 



The coexistence surface equation for n-player game is the 
Jacobian field for n-component output field, it is equiva- 
lent to the n-dimensional generalized Poisson bracket, 

..,*"}= ± ^§^...§^ = 0. 

(4.19) 

This coexistence equation of m vector field {<fii, (i = 
1,2, ...,m)} may be decomposed as a group equation of 
two field equations {{</>% <ft} — 0, i,j — 1, 2, m} where 
i and j must runs over all component of the vector field. 

In order to verify the universal coexistence curve 
equation, we take the two-phase coexistence equation 
{</> , <f> 2 } — as an example, and apply it to thermo- 
dynamic physics. The output field is the difference of 
free energy O u t("f) = SF = F A — F s , the input vector 
are temperature 71 = T and pressure 72 = P. It will 
be shown, the universal coexist equation {cp 1 ,^ 2 } = 
unified all the coexistence equations in classical phase 
transitions. 

We first verify the second order phase transition. The 
order parameter of the second order phase transition is 
6 1 = drSF and <fi 2 = dpSF, substituting them into the 
Jacobian vector 



{<j>\<l) 2 }=D{<t>lq) = 
and using the relations 



def) 1 d(j) 2 
dT dP 



dej) 1 d(j) 2 
dP dT 



0, (4.20) 



drdrSF 



- C, 



T 



d P d T 6F = V(a l 



B 

p -, dpdpSF = V(r4- k|), 
a A ), (4.21) 



we arrive 



V 



D(<p/q) = -(Cf - C A )(4 - 4) ~ (Va 



B 



Va A ) 2 . 
(4.22) 
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FIG. 3: The saddle surface of the free energy around the 
critical point. 



Recalling the Ehrenfest equations 



dP 
dT 



dP 



-Ct 



dT TV(a B ~ a A ) ' 



(4.23) 



it is easy to verify that the equation above is in consistent 
with the bifurcation condition Eq. (14.22[) . So the bifur- 
cation equation D(<f>/q) — is an equivalent expression 
of the coexistence curve equation. The solution of this 
equation is a two dimensional coexistence surface, the 
Ehrenfest equations actually indicates the normal vector 
of phase A and B is of equivalent value with opposite 
direction, i.e., \ub\ = — |^a|i that means they reach bal- 
ance on this surface. 

For the first order phase transition, we chose the vector 
order parameter as <j> = d°SF, here '0' means no deriva- 
tive of the free energy. The generalized Jacobian vector 
of the first order phase transition with 4> = d°SF is given 
by 



dC-) = { 



Q F B 

~df~ 



dF A ,dF B 
-) + (" 



dT 



dP 



dF A 
~~dP 



)=0, (4.24) 



in mind of the relation 



3F 



considering D(&) 



dT 

0, we have 



-S and 



dP 
dT 



(S B - S A ) 
(yB _ v a^ ■ 



V , and 



(4.25) 



This is the famous Clapeyron equation. The critical point 
is a saddle point. So it is the maximal point of free energy 
difference in 71 direction, and it is minimal point in 72 
direction, their first order derivative must obey 



dSF d5F 



< 0. 



(4.26) 



The first order phase transition requires they must share 
the same absolute value at the critical point, then 



d5F d6F 



dyi d^ 2 



0. 



(4.27) 



In fact, the free energy difference between the two sides 
of the coexistence acts as the phase potential, its first 
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derivative is force, the force of the two phases must be of 
the same value but pointing in the opposite direction at 
the critical point. 

The bifurcation equation D(<p/q) = can be natu- 
rally generalized to a higher-order transition, it also leads 
to the coexistence curve of the higher order transition. 
We consider a system whose free energy is a function of 
temperature T and magnetic field B, then the Clausius- 
Clapeyron equation becomes dB/dT = — AS /AM. If 
the entropy and the magnetization are continuous across 
the phase boundary, the transition is of higher order. For 
the pth order phase transition, the vector field is chosen 
as the (p — l)th derivative of SF, (ft 1 = d^~ 1 5F, <ft 2 — 
d%" 1 SF. Substituting (</>\0 2 ) into Eq. (|4~2T)|) . we arrive 



D(</>/q) 



qpSF d p 5F 



o-l OT^RP-I ^' ' ; 



DTP dBP dBdTP' 1 dTdBP- 



Considering the heat capacity 1^ = —^f- and the sus- 
ceptibility ^2 — Xj the bifurcation condition D(tfi/q) = 
is rewritten as 



dB 
dT 



(-1)1 



AdP- 2 C/dTP-'' 
T c AdP- 2 X /dBP- 



(4.29) 



This equation is in perfect agreement with the equations 
inRef. 

In mind of our holographic definition of phase transi- 
tion Eq. (|3.36p , we may also derive the holographic coex- 
istence equation using the fundamental order parameter 
field, 

2 S' P [R(9)(F A — F B )] _ 6P[R(0)SF] 



8 Op 



S6p 



R(9)SF 



SF, 



t d d 

L = ^i^— ■ 

Oil Ol2 



(4.30) 



To study the pth order phase transition, one need to ex- 
pand the group element R(9) to the p the order, and split 
it into the real part and imaginary part, i.e., 



6 p [R(6)6F] 
Wp 



9=0 



b^im 2 . 



Then the coexistence curve equation is 

I At A?\ ^ d ^ d(t)1 d ^ n 
071 0J2 Ol2 C71 



(4.31) 



(4.32) 



Under this definition, it is easy to see that Kunmar's re- 
sult (|4.29|) is only special case of a series of coexistence 
equations for the pth order of phase transition, the com- 
plete coexistence curve equations are given by 



d l 8Fd 2 P- l SF did p - j 6F d k d p - k SF 



g T i 9B 2 p- 
(i,j,k= 1,2, 



dBWTP-i 3T k dBP- k 



,P). 



0. 

(4.33) 
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FIG. 4: In analogy with Newtonian mechanics, the difference 
of output between two players is equivalent to gravitational 
potential V, the first order phase transition means the two 
phases have the same potential. For p = 2, the order parame- 
ter field, 4> A = 9^~ 1 O ut , 4> B = d^^Out, represents the force 
of phase A and B. When they reach balance, the two phase 
coexist on the hypersurface. This is the first order phase 
transition. For p — 3, <f) A and (f> B represents the accelera- 
tion, although the force is not equal, but the acceleration are 
the same, this is the second order phase transition. On the 
two sides of the coexistence surface, the tangent vector are 
continuous, but the normal unit vector has a sudden jump. 



Now we see, the universal coexistence equation not 
only reproduced all the coexistence equations of classical 
phase transition in physics, but also gave more equations 
that have never been appeared before. This indicates 
that the game theory of renormalization group transfor- 
mation has very broad applications. 

For the magnetic field and temperature depended free 
energy F(B,T), the scaling laws were derived in Ref. 
[l6l |. the exponents is defined as 



8p- 2 C 
OTP- 2 



qp- 



'X. 



dBP- 



(4.34) 



The equivalent expression in terms of the vector field 
is 



drd) 1 = 



dpcf) 2 = a 



(4.35) 



This scaling law is only a special case, there were many 
other scaling laws in various physical system, the same 
value of critical exponents falls into the same universality 
class. In the next section, we shall discuss the scaling laws 
based on the most general definition of phase transition. 



4.3. Universal scaling laws and coexistence 
equation around Nash equilibrium point 

The phase transition of a physical system occurs at 
the critical point where the correlation length between 



17 



particles becomes infinite [2]. It is assumed that the 
free energy is a generalized homogeneous function, i.e., 
F(A a 7 1 ,A f, 7 2 ) = AF( 7 1 ,7 2 ). The quantity defined by 
the free energy obey power laws around the critical point. 
From the two scaling exponent a and 6, one may derive 
those critical exponents which obey some equalities. The 
systems with the same scaling law fall into an universality 
class. 

When it comes to the general phase transition defined 
as a war game in this paper, all the critical phenomenons 
in physical reappeared. During the war, the correlation 
between the all the members of the participants becomes 
infinity, people may not know each other, but every tiny 
work they do may cause great effect to the final results. 
If we focus on the individual person in battle field, one 
would see two opposite soldiers are fighting. Then we 
go to a larger scale, we see two companies are fighting. 
We can continue to magnify the scale, the participants 
who are fighting range from hundreds of people to mil- 
lions of people, range from a small village to the whole 
world. No matter from which scale we see it, it is the 
same war. At the critical point, the participants of the 
war have equal strength, if any one of them make a tiny 
mistake (the mistake may come from an unimportant sol- 
dier), the whole army will lose the war, so the correlation 
length between soldiers goes to infinity. In this sense, the 
output field of the war should be a generalized homoge- 
neous function at the critical point, it obeys the relation 
0«t(A 4 7l ) = AO ut ( 7i ). 

According to the topological current of phase transi- 
tion, we established the tangent vector field of O ut . The 
tangent vector field is the projection of physical field con- 
figuration on the strategy manifold. The physical field di- 
vergent at some singular points where the tangent vector 
field vanished. It was proved that the nontrivial Rieman- 
nian curvature just around these surviving strategy. The 
integral of the Riemannian curvature is a topological in- 
variant, the critical exponent should bear a topological 
origin. 

According to the universal definition of phase transi- 
tion, the phase transition point is a Nash equilibrium 
solution. A special two dimensional output manifold is 
a saddle surface in the vicinity of the critical point. The 
output manifold is maximum for one parameter, but min- 
imum for another parameter. 

Suppose the topological dimension of the manifold 
around the critical point is integer, i.e., D — 1,2, ...,n, 
we may introduce a local coordinates to approximately 
express the output manifold in the vicinity of the critical 
point as 

Out = f 7? + f 7l, (4-36) 

where 71 is the infinite small variable, based on which the 
local coordination is r\ = (Sri + ^}j'^i'yj)y/gii + • ■ ■ • Wc 
may abandon the quadratic term of 7 j . Eq. (|4.36j) is the 
local approximation of the output manifold. k,\ and K2 
are principal curvature. When n\ > and K2 > 0, it is a 



elliptic surface, when k\ > and K2 = 0, it is parabolic, 
for Ki > and Ki < 0, it is hyperbolic. Usually the 
local manifold on a two dimensional output manifold is 
hyperbolic at a phase transition point. In fact, if we carry 
out the derivative of Eq. (|4.36|) to the seconde order, it 
spontaneously leads to, 

a?A t = Ki, d^6 ut = K 2 , (4.37) 

they are actually intrinsic geometric constant of the 
neighbor manifold around the critical point. 

Recall the game theory of renormalization group trans- 
formation in the first section, one would see that the 
game operator is a nonlinear operator, it defines an iter- 
ative map for the game process. The dimension of the in- 
finite small neighboring manifold around the Nash equi- 
librium point can be exactly calculated from the game 
operator. For most nonlinear game operators, the di- 
mension of the manifold around the Nash equilibrium is 
fractal instead of integer. There are only a few very sim- 
ple cases that one can find an integer dimension. But the 
game operator in that cases is too trivial to give us any in- 
teresting phenomena. According to the experiments and 
numerical calculation of physics, we can make a general 
hypothesis that the neighboring output manifold around 
the critical point of phase transition, namely around the 
Nash equilibrium point of a non-cooperative game, has 
fractal dimension. 

In the vicinity of the phase transition point, an approx- 
imation of the scale invariant output manifold is a com- 
plex function in fractal space, 

Outil) = Kitf + «m[/(7? + < k -)\ dm +■■■, (4-38) 

where c?i,d2, ...,d n are fractal dimensions wit respect to 
different input parameters. Recall that most of the phys- 
ical observables in statistical mechanics are defined by 
the second order derivative of free energy, we can define 
the observables of a complex system by the second order 
derivative of the output field. The free energy is only 
a special component of the output manifold in physics. 
When we study the most general complex system, as long 
as people can measure it, we can take any order of deriva- 
tive of the output field as observables. These observables 
are just the components of the tangent vector field of the 
output field. A simple example of the tangent vector field 
is 

^ = ~d^~ = daC + '"- (4 - 39) 

According to the topological phase transition theory, the 
tangent vector field satisfy the phase coexistence equa- 
tion {(f) 1 , <fi 2 , eff 1 } — at the critical point, we can de- 
rive a constrain on these vector field. We substitute the 
explicit expansion of the observable quantities into the 
coexistence equation, it would lead us to a constrain on 
the fractal exponent index. These constrain relations are 
just the scaling laws. Therefore scaling laws come from 
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the coexistence equation of the crossing defined physical 
quantities. 

The scaling relations found in various physical system 
are probably the simplest relations on the fractal space 
extended by two parameters. One can reach all kinds of 
different scaling relations in statistical mechanics by 
taking the output field as free energy O ut {ji) = F(7i), 
and taking the input ji as physical parameters, such 
as temperature T, pressure P, magnetic field B, and 
so on. In the vicinity of a second order phase tran- 
sition, the divergent physical quantity are defined by 
the second order derivative of the free energy. Such 
as the susceptibility \ — —d 2 F/dH 2 , H is magnetic 
field. Each divergent quantity is characterized by a crit- 
ical exponent, this critical exponent comes from fractal 
space. The two component coexistence {(p 1 , </> 2 } = pro- 
duced the scaling relations, such as the Fisher relation 
vd = 2 — a, Widom relation 7 = (3(6— 1), Rushbrooke re- 
lation a + 2/3+7 = 2, and so on. We take the Rushbrooke 
relation as example to verify the coexistence equation. 
The Gibbs free energy is G = U — TS, its differentia- 
tion is dG = -SdT + VdP - MdH. Experimental and 
numerical calculation found that three thermodynamic 
quantity obey the following scaling laws in the vicinity 
of critical point, 

M = -(§) ~ |Tf , 



Cr 



-T( 



d 2 G, 



\T\ 



Q T 2 

x = -(g) ~ ITT. 

The fundamental vector field can be taken as 



* = (—), 

K dH h 



6 2 = (— ). 



(4.40) 



(4.41) 



Substituting the two vector field into the coexistence 
equation 2 } = 0, one may derive 



d 2 G d 2 G d 2 G d 2 G 



dT 2 dH 2 dTdH dHdT 



0. 



(4.42) 



Now we substitute the thermodynamic quantities (|4.40|) 
into the coexistence equation, it yields 



ITI 



= /3 2 |T| 



2/3-2 



(4.43) 



When T — > 0, we may ignore the coefficient /3 2 at the 
right hand side of Eq. (j4.43[) . then we obtained the Rush- 
brooke relation a + 2/3 + 7 — 2- Other scaling relation can 
be verified following similar procedure. These relations 
were firstly found by computational simulation and ex- 
periments. Therefore the scaling law of universal phase 
transition in a general complex system has solid numeri- 
cal and experimental foundation. Here it must be pointed 
out that the commutable relation BtOh = 9h9t have 
been used in the calculation. This suggests that the par- 
tial differential corresponding to different variables are 



commutable in the vicinity of critical point. This is in 
consistent with our picture of war game at the phase 
transition point. Further more, one may choose different 
tangent vector field for the coexistence equation, then one 
may obtain all different scaling relations in the vicinity 
of critical point. 



4.4. The symmetry of Landau phase transition 
theory and symmetry of game theory 

The Landau theory of continuous phase transition the- 
ory provides a basic description to the phase transition 
characterized by spontaneous symmetry breaking. Take 
its application to the structure phase transitions as one 
example, it derives several important features, namely, 
the change in the crystal's space group, the dimension 
and symmetry properties of the transition's order param- 
eter, and the form of the free energy expansion. It has 
the same range of validity as the mean-field approxima- 
tion in microscopic theories. A central assumption of the 
Landau theory is that the free energy can be expanded 
as a Taylor series with respect to the order parameter 77: 

F(P, T, 7 1 ) = F + A(P, T) V 2 + B(P, T)rf + C(P, T) V 4 + ... 

(4.44) 

in which the phases are marked by the order parameter 77. 
The symmetry in Laudau theory talks about the invari- 
ance of this free energy when we do some transformation 
on the order parameter field 77' — > Urj. This symmetry 
is the same conception as that in quantum field theory, 
such as the tjj 4 model, C = ^(d^tp) 2 — im 2 i// 2 — jtp 4 . If 
the equations of motion derived from this Lagrange is in- 
variant under some transformation tp — > tp' , we call such 
a transformation a symmetry transformation. 

But in this paper, the symmetry we are talking about 
in the game theory of phase transition is a different con- 
ception. 

The object of our research is a very general system 
O u t = F ou t(x,j), 7 is the input vector. In game theory, 
the input vector represent the input strategies of different 
players. In physics, the input vector are physical oper- 
ational quantity. O u t is output vector, it encompasses 
all the information the observer can received by send- 
ing different inputs. In physics, the outputs are physical 
observables. The state x represent the inner states of 
system, it plays a similar role as the order parameter 
field 77 in the free energy expansion equation (|4.44ll . In 
our game theory of topological phase transition theory, 
all the state vector have been integrated out, the funda- 
mental starting point is the output field O u t = F out (^), 
the symmetry we mentioned in this theory is about the 
transformation invariant property of O u t under the trans- 
formation 7* — > U"f. For example, in the free energy 
equation (|4.44p . F(P, T) is output, P and T are the two 
players, we study the symmetry of F(P, T) when T — * T' 
and P — ► P' . The order parameter r\ is not an operation 
quantity, it describes the inner state of the system. 
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The Landau theory of phase transition can be sum- 
marized as a differential game(see Appendix IDIE[) . The 
order parameter is the state vector of the game. In the 
frame work of our topological current theory of phase 
transition, the tangent vector field <j) 1 consists of the com- 
plete set of phase dynamics system. The vector field 



— u p 



.! 6F(P,T,r,) 
( 5r) ' 



(4.45) 



describes how the two player T and P behave across the 
state vector space 77. The free energy is the output func- 
tion. F(P, T, ij) should be written as an expansion in the 
even powers in the spirit of Gintzburg-Landau formal- 
ism, for F(P, T, 77) must be gauge invariant with respect 
to the order parameter 77. Then the series of free energy 
in even powers truncated at the fourth order is 

F(P, T, r]) = F Q + A(P, T) V 2 + C{P, T)rf, (4.46) 

A special phase vector field may be chosen as 

1 = 2d^[A(P,T)+2C(P,T)r ) 2 }r ) , 

<f> 2 = 2dp [A(P, T) + 2C(P, T)r] 2 ]r]. (4.47) 

If (j) 1 and (j) 1 do not commute with each other, both the 
two players have surviving strategy which are the solu- 
tions of (j) 1 =0, <j) 2 — 0, these solution sit at some iso- 
lated points 



T k =T k {r{), P k = P k (r}), k = 1,2, 



(4.48) 



Each of these isolated solutions has a winding number 
W k . These surviving strategy pair (T,P) are changing 
according to the state vector rj. The surviving strategy 
of new phase and old phase carry opposite winding num- 
ber, they are annihilating and generating at the phase 
coexistence region. The sum of these winding number is 
a topological quantity which is determined by the topol- 
ogy of the state vector manifold. . 

Usually A(T,P) has the form of a(P)(T - T c ) near 
the critical temperature T Cl and C(P,T) is supposed to 
be weakly dependent on the temperature, i.e., drC <C 1. 
Considering Eq. (|4.47p , the surviving strategy for the old 
phase and new phase can be derived from 



r,[a(P) 
r,[d P a{P){T-T c ) 



2d T C(P,T)r) 2 }=0, 
- 2d P C(P, T)r] 2 ] = 0. 



(4.49) 



One sees that rj = corresponds to a stable phase solu- 
tion of the equation above. For the case 77 ^ 0, there are 
two solutions: 



77 = ± 



-a(P) 
2d T C 



1/2 



V = ± 



d P a(P)(T c - T) 



2d P C 



1/2 



(4.50) 

if we chose C = exp(T + P) and a = exp(P), it turns 
into a familiar form 



m = ± 



-a(P) 



2C 



1/2 



! m = ± 



a(P)(T c ~T)_y /2 

(4.51) 



2C 



The coexistence curve equation {(j) 1 , (j) 2 } = is an univer- 
sal equation, it holds for this two player game. Consider 
the two special vector field Eq. (|4.47|) . we arrived a so- 
phisticate coexistence equation, 

Arj 2 dlC[d 2 P A(T,P) + 2d 2 P C7] 2 } = [d P a{P) + 2dpd T CT 1 2 ] 2 . 

If we have obtained the explicit relation of A(T, P) and 
C(T, P), then we can depict the phase diagram in T — P 
plane for different state vector 77. 



4.5. Symmetry and evolution of phases in game 
theory 

In game theory, the output vector are payoff functions, 
the inputs are strategies. A symmetry transformation in 
strategy space does not always lead to different output, 
the Nash equilibrium solution is invariant fixed point un- 
der continuous transformation. As shown in Ising model, 
the renormalization group transformation keep mapping 
a pair of strategy to another pair, no matter where it 
starts, it always flows to the fixed point. 

The phase of a system evolutes in the direction of 
renormalization group transformation flow. The direc- 
tion of renormalization group transformation flow is de- 
termined by the Second Law of Thermodynamics, which 
states all physical systems in thermal equilibrium can be 
characterized by a quantity called entropy, this entropy 
cannot decrease in any process in which the system re- 
mains adiabatically isolated. We can ignore the thermo- 
dynamics, and grab the central point that an isolated sys- 
tem only evolutes in the direction of increasing entropy. 
This statement of the Second Law of Thermodynamics 
may also hold in game theory, but notice that the entropy 
here in our theory is not the conventional conception de- 
fined in statistical physics, since we are not manipulating 
the physical particles, it is the entropy of players 7 in- 
stead. In physics, the players are those physical observ- 
ablcs. An isolated system evolves spontaneously toward 
a maximal symmetry. 

The entropy of the multi-player game measures indis- 
tinguishability of players. Let's look at an ancient battle- 
field with two armies ready to fight, before the beginning 
of the combat, the two troops are organized in ordered 
states separated by vacant glacis, every part of them has 
special functioning. A bystander can tell which is which. 
Once the combat begins, they rush at each other and 
fused into one. In this case, the bystander is unable to 
distinguish them, it is completely in chaos. When the war 
finished, the loser surrendered, winner gathered, the two 
troops translated into another ordered states. The sur- 
viving troops are distinguishable again, but those dead 
can never be back to ordered states again. 

High symmetry leads to high entropy. If we do some 
transformation to a system in chaos, it won't make any 
difference. As in the war, the two armies become a mix- 
ture, it is an unstable equilibrium state with the highest 
symmetry. The higher symmetry a system owns, the 
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more unstable it would be. The state before and after 
the war are both of less symmetry states. Like the be- 
ginning of universe, all different conflicting interactions 
confined in a singular point, reach an unstable equilib- 
rium state. It has the highest space time symmetry. A 
minor imbalance results in the Big Bang. The explosion 
breaks the singularity apart into all different symmetry 
zone. Each symmetry group brings about a stable phase. 

The stable phases at different stage of history may be 
found following the group chain of the original symmetry 
group. A high symmetry group usually has a series of 
subgroup. A physical system liken to stay in the most 
fundamental symmetry state. People love peace, but hate 
war. They like to live at the most stable point of the 
system. A war cost too much energy. The subgroup 
chain of a given group is always finite, 



U DU!DU 2 D U 3 ... D U p . 



(4.52) 



For example, 5Q(4, 2) constitutes a dynamical group for 
the Hydrogen atom[l5}. It has many four subgroups [l7| . 

5*0(4,2) D 50(4, 1) D SO(4) D 50(3) D 50(2). 

(4-53) 

This group chain may be followed by point group. For 
instance, the U(l) symmetry group may be decreased 
to T n which means a physical system is invariant when 
it rotates ^ around one axel. So we can roughly find 
the stable phase according to the group chain. Generally 
speaking, Lie group is highes rank of group. The discrete 
group are imbedded in Lie groups. The most stable phase 
focus on the lowest symmetry. 

A game system evolutes following the first principal, 
the players take different strategies to ground state. The 
N players represents N interaction, or N input vectors. 
In the beginning, each of them occupies a small domain 
in phase diagram and is the ruler in his own domain. 
There is temporary balance between the different play- 
ers. The war is going on all the time between neighboring 
domains. Whenever a player loses, he dies out, his ene- 
mies absorbed his domain. When a player fails, an old 
phase died, and a new phase is born, this indicates a 
phase transition. The balance on the boarder between 
different domains represent the phase coexistence region, 
on which Nash equilibrium is reached. But this state is 
an unstable equilibrium state, there is constantly minor 
conflicting to destroy this balance. Then the war goes 
on, the winners gain advantage, and becomes stronger 
and stronger, the war will not stop until he unified the 
domains as one. 

The phase diagram of a physical system is depicted 
in the frame of some static physical parameter. We can 
view the evolution of phase as an continuous tuning of 
the physical parameter. For example, if we continuously 
raise temperature, we will see single curve in the phase 
diagram bifurcates at certain temperature, this is the 
evolution of physical phase, although it is not the usual 
evolution in the common sense of time. 

This phase evolution with respect to temperature can 
be viewed as a generalized dynamic process in parameter 



space. As shown in section (|3.2[) . a multi-player game can 
be equivalently expressed as a group of differential equa- 
tions. The evolution of output vector O ut is governed by 
the differential equations, 



d s O ut = X(7,<9 7 )O ut , i = 1,2, . 



(4.54) 



If we take the tuning parameter s as temperature T, then 
we understand the evolution of phase with respect to 
temperature. In physical system, we can obtain the func- 
tion of linear response by perturbation theory. The re- 
sponse function O u t may have sophisticate relation with 
input parameters, O u t(^). Some input parameters are 
function of temperature ji = 7t(7j). For example, dur- 
ing the isothermal process, the product of pressure and 
volume is a constant PV — constant, so P and V are 
related by this equation. We take the as funda- 

mental variables, jj is the tuning parameter. Then the 
generalized momentum is defined as 



P, 



0(a»r) • 



(4.55) 



where 4>('J, <9 7j 7, 7j) is a reference function like La- 
grangian function. We can further find a vector field 
4>h{j, 9 7j 7, 7j-) which is similar to the Hamiltonian func- 
tion in Lagrangian mechanics. For a given field 4>h, the 
game operator £(7, <9 7 ) bear an expression in terms of 
Poisson bracket, 



L(7, d 1 



d 



d 



dP 7i d"f c?7 <9P 7i 



(4.56) 



This formula of game operator leads to another equiva- 
lent representation of the differential equation , 



d lj O ut = {4>H,O ut }. 



(4.57) 



This equation actually describes the game process be- 
tween <f>H and Out, each of them navigates one stable 
phase. If they are in the same complete set, O u t} = 
0, then the two phases coexist. In fact, the Hamiltonian 
function does not play any special role in our topological 
current theory of phase transition, it is just one of the 
tangent vector field in Lie algebra space. 



4.6. Phase coexistence boundary and unstable 
vacuum state 

As we defined in the beginning, phase transition is a 
transition from one stable state to another stable state, 
there is a critical point at which the old stable state col- 
lapsed and the new state arise from the shambles and 
grow to a stable state. A oversimplified model is take two 
states as a vector of two component (2, y), the evolution 
operator is a diagonal 2 by 2 matrix L = diag(— 1,1) 
at the coexistence region. The evolution of the two 
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phase follows the equation <9 7 r = Lr, which we estab- 
lished in last section. Then in the vicinity of the coexis- 
tence region, we would see the old phase decays follow- 
ing x oc e _ ° 7 , in the meantime the new phase blows up 
y ex e +Q7 , where a > 0. 

The mechanical simulation of the phase coexisting 
state is a ball on the maximal tip of the parabola f(z) = 
— z 2 . At Nash equilibrium point, a player take the best 
strategy to minimize his own damage, and obtain his 
maximal profit in the meantime. In fact, it is the profit 
difference between them that they fight for. When the 
two conflicting force reach a balance, the game arrived 
at a Nash equilibrium. 

The Nash equilibrium state is an unstable maximal 
point in the mechanical potential of the output field. Let 
% u t{l) be the output function of the player 7*, the ef- 
fective potential for two players in a game is AOy = 
Outil) ~ Outil)- Nash equilibrium sits at the minimal 
point of |<50y|. Notice that the player 7"s profit is the 
damage of player 7-?, so AOij — — AOy. We may es- 
tablish a general equation of motion for player 7*'s profit 
function, 

d 2 y & ut =d r V(& ut ), (4.58) 

V(O l ut ) is a general potential which is self-consistently 
decided by other players, usually the player 7* sits at 
some of its minimal points, while the other players occu- 
pied its maximal point. The typical profit function O l ut 
of player 7* is a kink solution, 

6^±tanmy-7*)= e lil .T e l-1 - ( 4 -59) 

e 7 7 4- e 7 7 

where 7* is the optimal strategy. Eq. (|4.59|) is the fa- 
miliar kink solution of quantum tunnelling problems fl8l|. 
Any minor deviation from the Nash equilibrium solution 
would results in drastic increase or decrease of the profit 
function. Therefore the Nash equilibrium solutions of a 
game play the same role as the vacuum solution to quan- 
tum tunnelling. A collection degenerate vacua means 
there are a series of Nash equilibrium solutions with the 
same optimal value. 

The game of phase transition has two different equi- 
librium states. One is all different interactions reach an 
agreement to maintain peace, this is the optimal strat- 
egy of cooperative game. It is the trivial vacuum. The 
other is that no peace agreement is derived, a war breaks 
out. At the critical point, different phases coexisted but 
against each other. This is the Nash equilibrium state. 
It is an unstable vacuum state. 

The coexistence equation of n-different phases we ob- 
tained in previous sections, 

{0\0 2 ,...,<n = O, (4.60) 

describes a Nash equilibrium state, or in other words, 
the coexisting unstable vacuum state, here <f> = 



(0\ 2 , (j) n ) = S " [U ^l 0ut] \e=o. For this equation is an- 
tisymmetric, exchanging any two of the players would 
add a (—1) to the output. As shown in the topologi- 
cal current of phase transition, the sum of the winding 
numbers around the surviving strategies is a topological 
number — -Chcrn number. The sign of each winding num- 
ber is determined by the sign of {(ft 1 , (f> 2 , 4> n }. Every (jf 
field is a player, if exchange any two of them, the winner 
becomes loser, the loser turns into winner. Therefore if 
(j) 1 has a positive winding number, his opponent must has 
a negative one. The strategy for (j? to survive is the anti- 
strategy of his opponents. During phase transition, the 
old phase is the opponent of the new phase. The wind- 
ing number of the new phase's strategy 7 J is Wjf ew > 0, 
the winding number around the surviving strategy j l of 
the old phase is W % old < 0. The topological constrain 
suggests 

NchernNumber = ^ W Old + ^ W NeW ( 4 ' 61 ) 
i 3 

We may view each strategy of the new phase as one parti- 
cle with topological charge Wn £W , and its corresponding 
anti-particle is the anti-strategy of the old phase which 
carries a negative winding number. The soldiers of the 
new phase are the strategies, they carry positive wind- 
ing number, rush at the coexist curve to fight against the 
anti-strategy of the old phase. The particles and antipar- 
ticlcs annihilate at the phase boundary, so the coexisting 
phase boundary behaves as vacuum. This war is going 
on under the constrain of Eq. (|4.61|) . 

The generator of translation group along 7^ is id 1 ^. 
The evolution of the output on the strategy space follows 
the classical Hamiltonian-Jacobi equation, 

id^{6 ut {9)} ={{O ut {9)), L}, (4.62) 

The second quantization of this equation is just the 
Heisenberg equation id t O u t = [O u t, H], which has a more 
general covariant form 

id^6 ut = [d uU P^ (4.63) 

here 7^ = (£,7), P^ = (H,P). In four dimensional 
Minkoveski space time, the Hamiltonian is the generator 
of translation group with respect to time t. While the 
momentum operator is the generator of space trans- 
lation group. Integrating the covariant Heisenberg equa- 
tion (|4.63p . we arrive 

0( 7 ) = e J ^'O(0)e- l/>7 ". (4.64) 

If the output is expressed by quantum operators, its evo- 
lution in strategy space is governed by renormalization 
group transformation 
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U(d){0)U-\6) =e Wn - h {0)e 



-iOfi-L 



(6)+i6 [n-L, (6)} + 
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{i6f 
2! 



[n ■ L, [n ■ L, (O)}} + 



(4.65) 



For the most simple Lie group 50(2) whose Lie algebra 
has only one generator L z , this transformation equation 
takes a very simple form, 



us(6)u- 1 

= 6(6) + 18 [L z , 3(6)} + ^-[L z , [L z , 5(6)}} + 



[i6f 
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In quantum mechanics, if a group transformations U 
commutes with Hamiltonian H, [U,H] = HU — UH = 0, 
i.e., H = H' = UHU\ then H is invariant under trans- 
formation of group U, they share the same eigenfunction. 
In a general game theory, the Hamiltonian has nothing 
special. The generators of Lie algebra plays the role of 
Hamiltonian operator. The quantization of the coexis- 
tence equation reads 



0. 



(4.66) 



Each operator <p l represents a quantum operator, if 
they commute with each other, they share the same 
eigenspace. Eq. (|4.66|) is the quantum coexistence equa- 
tion. 

If the 4> % are vectors expanded in the tangent space 
of Lie group around the identity, they are vectors of Lie 
algebra. For Mathematician, this coexistence equation 
corresponds to an invariant Cartan space of Lie alge- 
bra. The coexistence phase space is the eigenspace of 
Cartan subalgebra. The Cartan subalgebra is the maxi- 
mal Abelian subalgebra. An arbitrary Lie algebra vector 
commuting with this subalgebra is still a Lie vector in 
the same space. The elements Cj in this commutative 
subalgebra must satisfies 



[c i ,c j } = o, (ij = o,i, 2, ...,/). 



(4.67) 



I is the dimension of the Cartan subalgebra space. It is 
also the number of coexistence phases. For the coexist- 
ing point of three phases, we define the three operator 
commutator as 

[Cj, Cj, Ck] 
= [Ci,Cj}Ck + [Cj,Ck}Ci + [Ck,Ci}Cj 
= 0, (i,j,k = 0,l,2,...,l). (4.68) 

It has a straight forward generalization for the points at 
which n phases intersects, 

[Ci,Cj,C k , C f ] = 0, k, / = 0, 1, 2, ...).(4.69) 

So we have a well defined quantum operator 

D" = [^ 2 ,...,0 n ], (4.70) 



This operator is defined from by vectors of Lie algebra, 
we call it phase coexistence operator. According to the 
group representation theory [20j, we can always find the 
representation of coexistence operator, 



D n \^ n ) =£>|V„). 



(4.71) 



If n = 2, D n is just the commutator of two operators 
which represent two stable phases. The phase boundary 
between phase A and phase B is given by the zero modes 
of the commutator of the two phase operators. When 
(/>A and <Pb are not commutable, the eigenfunctions of 
the phase coexist operator may be divided into positive 
modes and negative modes, 



D n \il> n ) = signW\i) v 



(4.72) 



in which W is the winding number around the surviving 
strategy. sign(W) = +1 > corresponds to stable phase 
A and sign(W pt ) = — 1 < corresponds phase B, while 
(Wpt) = indicates the phase boundary. One of the most 
familiar example is to take the phase operator as the 
three component of angular momentum operator, cf>A = 
L x and 4> B = L y , then D 2 = L z . L z \ijj m ) = m\i/) m ), m = 
{0, ±}. m = is the coexistence eigenvalue, m = ±1 
represents the two stable phases. 

In topological quantum field theory, the topological in- 
formation of the zero mode is described by the Atiyah- 
Singer index theory. The strategy space of the game of 
phase transition can be split into positive eigenspace and 
negative eigenspace. The surviving strategies of the new 
phase carry positive winding number, we call them pos- 
itive modes </> + . The strategies of the old phase is the 
enemy of new phase, they carry negative winding num- 
ber, so we call them negative modes 4>~ . As shown by 
Eq. (|4.9[) . our topological current theory of phase transi- 
tion proved that the Euler characteristic number on two 
dimensional strategy space is, 



Ch = J <5^+){>+,<}- / S 2 (^){^,^}, 

= H w i-ll w -- ( 4 - 73 ) 

i 3 

For the 2n-player game, the topological Chern number 
on this 2n dimensional strategy manifold is determined 
by Atiyah-Singer theorem, 

Ch Bz = IndexD n = dimb n + - dimDI, (4.74) 

where D n is th quantized phase coexistence operator. 
The positive modes behaves as particles, and the negative 
modes are antiparticle, they coexist in vacuum. They are 
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born by pairs from vacuum, and annihilate by pairs to 
vacuum. 

In the language of Atiyah-Singer index theorem(see ap- 
pendix section U), the topological index is the difference 
between the number of positive eigen-modes and nega- 
tive eigen-modes. The positive modes is the surviving 
strategy for the new phase, the negative modes is the 
surviving strategy for the old phase. Thus the topologi- 
cal index counts how many extra strategies the new phase 
have after his soldiers rushed at front and annihilated at 
the coexistence phase boundary. We first assume every 
surviving strategy has a topological \W\ = 1 for con- 
venience, the Euler number Ch = +2 on a sphere, this 
topological constrain says that when all the soldier of the 
old phase died, there are at least two positive soldiers of 
the new phase left. The victory of the new phase has been 
determined by topology of the strategy manifold. If it is 
on a torus, the Euler number is zero, Ch = 0. The new 
phase is not that lucky now, his total number of soldiers 
is identical with that of the old phase. When the war 
break out, the strategy-anti-strategy pair annihilate on 
the phase coexistence boundary, the final result is a draw. 
When the base manifold of strategy space is a torus with 
n-holes , the Euler number is Ch = 2(1 — n), n > 1, the 
final victory goes to the old phase, he has at least 2(1 — n) 
surviving soldiers when the new phase is extinguished. 



5. PHASE TRANSITION AND 
ENTANGLEMENT IN GAME THEORY 

5.1. Phase transition and entanglement 

The participants of a war entangled with each other 
before the war breaks out, otherwise there will not be a 
war at all. In fact, a game is played by many interacting 
players. During the game, a player choose strategy ac- 
cording to other players' strategy, so he could never be a 
free particle unless he is not in the game. 

The entanglement we are talking about here is a much 
more general conception, it includes the relation, inter- 
action or connection between the elements of a system. 
The quantum entanglement in physics is one special case. 

When the players are at peace, the entanglement be- 
tween them is kept at a stable level. As the imbalance be- 
tween different players increase, the war is coming, their 
entanglement grows stronger and stronger. The whole 
system becomes more and more unstable. When the 
critical point arrived, a negligible event trigged the war, 
which spread the whole system through the strong en- 
tanglement. During the war, all the players summon up 
their internal strength to collect information from all the 
other players, and make the best strategy to win the war, 
the entanglement reach a climax. The external response 
during the war also reached the strongest level. When 
the war is over, everything is in order, their entangle- 
ment gradually decay to another stable level. 

Therefore, phase transition occurs when the entangle- 



ment between different phases reaches a maximal point. 
Every local maximal point of the entanglement indicates 
a transition, or a war. In order to give an exact pre- 
diction on where or when the war arise, we need to find 
some detectable quantity to measure the entanglement. 
So that we can quantitatively determine the position of 
the phase transition point. 

The most familiar quantity for physicist is the von Neu- 
mann entropy. Quantum statistics suggests that the von 
Neumann entropy of a pure ensemble is zero. Entropy is a 
quantity to measure how disorder a system is. More dis- 
order means higher entropy. The von Neumann entropy 
measure the disorder of the mixed states. The von Neu- 
mann entropy is a relatively small quantity if the system 
is in a stable phase. When the phase transition occurs, 
it would reach a climax point. 

In fact, any sensitive output function corresponding to 
a group of inputs can be used as measure of entangle- 
ment. If the inputs have strong relation between them, 
one minor change would definitely change the others, and 
this would leads to the response of many outputs. One 
can see this from a war, it is at the war that the ignor- 
able individuals began unite to work as a team, they are 
strongly correlated and entangled. 

In our topological current of phase transition, a quan- 
tity to measure entanglement may be defined as 
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where D((p/j) — {(f) 1 , <fi 2 , cjf 1 } is the coexistence equa- 
tion, and h is step size of renormalization group trans- 
formation. As shown in the previous section, the phase 
evolution equation in parameter space is d 1:j O u t = 

{<f>H,O u t} — LOut- So Eq. (|5.ip is some kind of prop- 
agator of phases similar to the propagator in physics 
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We rewrite the coexistence equation 



D{4>H) as Det where $ is the matrix extended by 
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Recall the definition of Pfaffian for matrix M, 



[Pf M] 2 =Det M, we can decompose the coexistence 
equation as D{(j)/^)=\Pi <f>] 2 . The entanglement equa- 
tion (|5.ip reads, 
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A typical physical example of Pfaffian is the fermionic 
parity. In the game theory here, a Pfaffian is a sum over 
all partition of the players into pairs, exchanging any two 
of them contributes a minus sign. For a bargain game 
between buyer and seller, the step size h is amplitude of 
the unfixed money they are fighting for. The smaller h 
they are negotiating on, the more information they are 
exchanging so that they could persuade each other to 
accept his offer. This means the entanglement between 
them increases when h approaches to zero. So the entan- 
glement increases when the renormalization group trans- 
formation goes on. The entropy increases in the mean 
time. The maximal entanglement appears at the phase 
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coexistence state, 



Eent = lim - !— 



exp | 
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<y(Pf$). (5.3) 



Here Pf $ = is an equivalent expression of the coexis- 
tence equation. The matrix $ in B e n< is constructed by 
output vector field. The output vector <j> could be any 
physical parameters. For the most simple case of two 
phases, <j> has two component <p A and <fi B , the Pfaffian is 
just the Poisson bracket. 

Besides the statistical observable or external response, 
the output vector can also be chosen as the conventional 
order parameter which is a quantity used to characterize 
the structural and inner order changes of physical system 
at the phase transition point. In the superconductor- 
insulator transition, upon tuning some parameter in the 
Hamiltonian, a dramatic change in the behavior of the 
electrons occurs, the order parameter of this quantum 
phase transition <f> = Ae l6 is the energy gap function of 
cooper pair theory. For the ultracold Bosonic atom gas 
confined in an optical lattice, the order parameter is the 
mean field value of the operator of Bosons <p = (b) = (¥). 
There were some numerical calculations in the Hubbard 
model, it shows the entanglement follows different scal- 
ing with the size on the two sides of the critical point 
denoting an incoherent quantum phase transition 21 1 . 

Entanglement state is the most important resource for 
quantum information technology. The entanglement be- 
tween the different stable phases is a good candidate for 
quantum computer. The most entangled states exist at 
the critical point of phase transition. The critical point is 
where the war breaks out, one of the fundamental char- 
acter of war is chaos, thus any minor change of parameter 
would leads to totally different results. Quantum com- 
putation using entangled states is extracting information 
from chaos, and control its output in a exact way. In 
fact, no matter it is classical system or quantum system, 
entanglement is the basic source of information manip- 
ulation, the best entanglement states exist in a chaos 
state. The biggest problem for present quantum comput- 
ing schemes is decoherence of entanglement, the quantum 
entanglement decays rapidly as time goes on. 

However, the entanglement we present in this section 
is independent of time, it only relies on the physical pa- 
rameter. It would be much easier to tune the physical 
parameters that to fight against time. So the entangle- 
ment in the vicinity of phase transition is very promising 
candidate for quantum computation. 



5.2. Quantum states and Nash equilibrium solution 
of games 

We have shown that phase transition occurs at 
the Nash-equilibrium point of a non-cooperative game. 
While the most entangled states just arise from this 
Nash-equilibrium point. We need to get a deeper under- 



standing on how the entanglement grows stronger and 
stronger in a game process. 

The most fundamental principal for a physical system 
is the Principle of Least Action. Newton's mechanics is 
unified in Hamilton's principle of least action as well as 
in Gauss's principle of least constraint. Maxwell's equa- 
tions can be derived as conditions of least action for elec- 
tromagnetic field propagation. When light goes through 
optical systems, it always take the path of least time, 
and takes short cuts in glass and water where light trav- 
els slower. The most stable state of many body system is 
the ground state, which possesses the least energy. The 
ground state is the final game results of the particles in 
many body system. Every particle wants to stay at the 
most stable point, it is the medal all the other particles 
struggle to get, this drives the particles in the game. 

The Principle of Least Action of a cooperative game 
is to maximize the profit of the whole group. The Prin- 
ciple of Least Action of a non-cooperative game is that 
all players only take strategies to maximize his payoff 
function or to minimize his loss function. Of course his 
personal strategy must take into account of other play- 
ers' strategy, because if the whole group breaks down, he 
will get nothing. 

The players may be classified by their statistics in anal- 
ogy with the statistics of particles in physical system. As 
all know, there are two types of elementary particles: 
Fermions and Bosons. There are anyons whose statistics 
stands between fermions and bosons, we first put that 
aside for convenience. The combination of fermions may 
form boson. Particles obeying Bose statistics shows a sta- 
tistical attractive interaction. While the fermions demon- 
strated a statistical repulsive interaction which comes 
from the Pauli exclusion principle. 

The predetermined physical environment is the game 
rule of the players. Players choose their states accord- 
ing to its interaction with external field and other play- 
ers. We can alternatively specify each strategy profile 
(si, S2, sjv) by the occupation number (rii, ri2, tin), 
which means there are rij players take strategy Si. Since 
all men are born free and equal, there would be no payoff 
difference if they choose the same strategy. 

The players in a non-cooperative game behaves as 
fermions in physical system. If we exchange any two of 
the players, the difference of profit between them would 
change a sign. They do not share the same occupation. 
So the expected number of players in noncooperative 
game is fermi-dirac distribution, 
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where Ui is the benefit the player could get by strategy 
state Si. The players in cooperative game are altruistic 
players, they concern for the welfare of others to maxi- 
mize the benefit of the whole group. They do not care 
too much about their own profit, so exchanging two of 
them does not make any difference, they like to share 
with other players. In this case, the distribution obey 
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Bose-Einstein statistics, 
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However, it is naive to say a player in reality is altru- 
istic player who only concerns about others or the ego- 
istical player who only concerns about himself. Every 
real player is in a mixed state of the altruistic and the 
selfish. We consider a rational player, each time he sets 
a step further by choosing a particular strategy, he has 
to confront a wining result or a losing results. When he 
loses, his benefit is transferred to other players, we call 
him altruistic, on the contrary, we say he is selfish. We 
expressed the altruistic state as |1) and the selfish state 
as |0), 



\altruistic) 
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(5.6) 



|1) and |0) form an orthogonal basis for the self-state 
space of player, i.e., (1|0) = 0. Because when he faces 
a particular pure strategy, there is only two possibilities: 
take it or give it up. An arbitrary mixed self-state vector 
is the linear combination of the two basis, 



|m) =a|l) + 6|0), 



(5.7) 



where a and b are complex number. |m) is a unit vector, 
(to | to) = 1, it is equivalent to a 2 + b 2 = 1. |to) means 
the player has the possibility of a 2 to be altruistic and b 2 
to be selfish. 

A game consists of N players is physically equivalent 
to a many body system with ./V particles. Each player 
represents a particle, the self state is the quantum states 
of the particle. A statistic distribution of the N particle 
states can be present in an ensemble which is a collection 
of identically prepared physical system. For a n player 
game, the global strategy state space is the product space 
of the n players strategy space, 



S = S" (1) <8> 5 (2) 
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The payoff function is a set value map from this Hilbert 
space to a number in Euclidean space. The payoff func- 
tion plays the role of a negative Hamiltonian in may body 
physical system. 

We first take the prisoner's dilemma to illustrate the 
basic phenomena of quantum states in non-cooperative 
game theory. 

In conventional text books about prisoner's dilemma, 
the two players are rational players(see Appendix [Gf . 
When they are prevented from cooperation, both of them 
would confess to minimize his own loss, both of them 
spend b years in jail. If Alice and Bob communicate and 
cooperate with each other, they would not confess so that 
they only serve a < b years in prison. 

Here we take a different angle to read the Prisoner's 
dilemma for the sake of quantum statistics. The self- 
states of two players have four different cases: (1) Al- 
ice and Bob are selfish players; (2) Alice and Bob are 



altruistic players; (3) Alice is selfish player and Bob is 
altruistic players; (3) Alice is altruistic player and Bob is 
selfish players. Selfish player take strategy to maximize 
his own benefit, and altruistic players tries to maximize 
the other's benefit. 

If we know which kind of players Alice and Bob are, 
we can find the Nash equilibrium fixed point, the two 
players does not need to commute with each other. For 
case (1), Alice knows Bob would certainly choose confess 
that can maximize his benefit, the only Nash equilibrium 
solution is that she also confess, so they reach the Nash 
equilibrium point (a, a). For case (2), the Nash equilib- 
rium is neither of them confess. In case (3), Alice does 
not confess, Bob confess. In case (4), Alice confess, Bob 
does not confess. The above is the case when the two 
players are in pure states, Alice(bob) is either altruistic 
1 1) or selfish |0). We denotes the four mixed self-states 
as 



\Da\D 



\1)a\0)b, \0)a\1)b, \0)a\0)i 



(5.9) 



The Nash equilibrium solution of the four states are sum- 
marized in table (|5. 10f> 
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(b, b) 
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This table is the density matrix of the prisoner dilemma. 
The statistical weight of |1)a|1)b is the Pareto optimal 
point of collective payoff. The statistical weight of self- 
state |0)a|0)b is Nash equilibrium. For the other two 
states |1)a|0)b and |0)a|1)b, the two equilibrium points 
appear in the off diagonal elements of the payoff matrix. 

The above is the simplest case of prisoner's dilemma, 
there are only two available strategies: confess or non- 
confess. Now we consider a more complex case: the avail- 
able strategies of Alice is {e^i, &ai, eAg} and Bob's 
strategies are {e^i, es2, &Bp\- The payoff matrix is 
now a q x p bi-matrix. 







Bob 




Bob 












Alice 




[An, B n ) 




(Aip, Bi p ) 












Alice 


&Aq 


(A q i, B q i) 




(Agp, B q p) 



(5.11) 



The prisoner dilemma told us, the fixed point solution 
of a game can be determined by self-state of players, 
namely, they are altruistic or selfish. A player could 
be in a mixed state of altruistic and selfish. For each 
pure strategy |e fc ), we introduce a fractional number pk 
to measure the altruistic degree of a player, pk satisfies 
the constrain J2k Pk = !• F° r example, we divided the 
altruistic degree into n stages, 



|0>, 



I) 



(5.12) 



26 




» 

o ° 

» 

FIG. 5: This is the payoff table which is rotated by 45 degrees 
in clock wise direction. Each black dot represents a pair of 
payoff value, the elliptic black point in the middle represent 
the Nash-Equilibrium(N.E.) points. The game begins at an 
arbitrary pair of strategies, there is an imbalance between 
Alice's payoff and Bob's payoff in the begging. Then the losing 
player takes a better strategy to eliminate the imbalance in 
the next step, the winning player also choose his best response 
to keep his predominance. This game process is actually a 
renormalization group transformation, it finally converges at 
the Nash-Equilibrium point. 



From |0) to |1) record how a player grows from a self- 
ish player to a altruistic player step by step. Once we 
know the self-state of a player, we know the probability 
he would play a certain strategy. As shown in Prisoner 
dilemma, the two player is either altruistic or selfish, the 
probability of confess as well as not-confess is either 
or 1. A vector of this altruistic degree could be the self- 
characteristic vector of a player, the equilibrium point of 
a game is utterly relies on this self-characteristic vector. 
There exists a one-to-one correspondence between the 
self-characteristic vector of players and the fixed point of 
a game. 

This self-characteristic vector of one player uniquely 
determined the probability he choose certain strategy. 
Thus when we introduce the mixed strategy for two play- 
ers Alice and Bob, 

<? 

k=l k 
P 

M = £pukl4>> 5Z/9 Bfc = l. (5.14) 

fe=l k 

the self-characteristic vector of Alice and Bob have 
been uniquely defined by {pAk} and {psk}- So every 
mixed state has an unique fixed point in the payoff 
matrix. In physics, payoff function or loss function 



corresponds to physical observables, such as Hamil- 
tonian, angular momentum, the given mixed strategy 
is just the cigenfunction with the fixed point solu- 
tion as the eigenvalue. We can define the density 
matrix using the strategy vector \s), p = Tr(\s)(s\). 
The von Neuman entropy S(p('j)) = —Tr(p(^) Ig p(-y)) 
measures the entanglement between the strategy vectors. 

Game theory of Hohenberg-Kohn theorem 

A many electron system can be mapped into a n 
player game, each player carries ±5 spin. The ground 
state energy of an electronic system is completely 
determined by the minimization of the total energy 
as a functional of the density function. The exter- 
nal potential together with the number of electrons 
completely determines the Hamiltonian, these two 
quantities determine all properties of the ground state. 
Hohenberg-Kohn theorem states: the external potential 
V is determined, within a trivial additive constant, by 
the electron density. This theorem insures the that there 
can not exist more than one external potential for any 
given density. 

The electrons are players, the many body wave func- 
tion are the strategy space. The external potential comes 
from external constrain. For example, in the prisoner 
dilemma, the police could put the two prisoners together 
or separate them. If the two prisoner are put together, 
they would conspire to remain silent so that both of spend 
less time in jail. If the police put them in separated 
rooms, they can not communicate with each other. Since 
they do not trust each other, they would confess. So we 
see external potential from police decided the behavior 
of the player, or in another point of view, the external 
potential transform the altruistic into the selfish by sepa- 
rating them, and transform the selfish into the altruistic 
by putting them together. Therefore there is only one 
self-characteristic vectors for a given external potential. 
This is Hohenberg-Kohn theorem in game theory. 



5.3. Renormalization group transformation and 
quantum entanglement 

In fact, the von Neuman entropy measures the entan- 
glement between the self-characteristics vectors. Every 
self-characteristics states vector is a mixed states of al- 
truistic and selfish. If the player is in a complex mixed 
states of altruistic and selfish, it is hard to predict where 
he ends. The Hohenberg-Kohn theorem told us the ex- 
ternal potential uniquely determined the state vector of 
players, as shown in prisoner dilemma. This suggests us 
a possible way to operate entanglement state using ex- 
ternal potential. 

The entanglement states grows stronger and stronger 
following the step of renormalization group transforma- 
tion. For example, in the bargain game, the seller first 
gave a price to see if the buyer take it. The buyer thought 
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it too expensive, he feedback his price to the seller. Then 
the buyer and seller both know each other at the first 
step. As this bargain goes on, they know each other bet- 
ter and better. In other words, the entanglement between 
them grows stronger and stronger. When this entangle- 
ment reach a maximal point, a Nash equilibrium state is 
arrived. At this point, if the seller raise one penney, the 
buyer won't buy it, on the other hand, if the buyer lower 
the price one penney, the seller will not sell his product. 
When the game passed over the Nash equilibrium, the 
entanglement decreased dramatically. 

The entanglement in non-cooperative game is much 
stronger than the entanglement in cooperative game. In 
last section, we know if the self-character of players is 
determined, we can uniquely find an equilibrium solution. 
We call the altruistic player an angle player, while the 
selfish a devil player. The angel player who manifests 
goodness, purity, and selflessness, behaves like bosons. 
While the devil player are fermions. 

The devil player do not trust each other, they most 
likely to betray in the game. So two devil players would 
try their best to bound them together and form a pair, 
they increase communication and cooperation to prevent 
his companion from betray. The entanglement between 
them is very strong, this entanglement would reach a 
maximal when the trial is around the corner. 

But angle players always trust their companions. Ev- 
erything they do is to increase the welfare of other play- 
ers. The angel player does not need to communicate too 
much. They behave much like non-interacting, indistin- 
guishable particles. So angle player is boson. An un- 
limited number of bosons may occupy the same state at 
the same time. At low temperatures, bosons can behave 
very differently than fermions; all the particles will tend 
to congregate together at the same lowest-energy state 
to warm each other, this is Bose-Einstein condensate. 

For a given many particle system, the entanglement 
can be measured by the particle-particle correlation 
length, the longer correlation length there exist be- 
tween particles, the stronger entanglement they have, 
we can read the entanglement out following Kadanoff 
block transformation. For example, in the two dimen- 
sional Ising model, the spins are players of game. They 
are mixed type of players between angle and devil. Their 
self-characteristic state vector can be expanded by 
four entangled eigenstates of the Bell operators: ■0 ± = 

^(| |)| t)±l T)l I)), ^ = ^(1 4)1 4)±| T)l T»- 

The density matrix, describing all the physical vari- 
ables accessible to entanglement state is given by 
p = Tr(\^p) (tpl). Then we can calculate von Neuman 
entropy S°(p('y)) — — Tr(p('j) Igp^)) for the first or- 
der Kadanoff block transformation. Then we construct 
the entanglement states between two blocks and obtained 
the second order von Neuman entropy S . We can con- 
tinue this renormalization group transformation, and fi- 
nally derive an exact von Neuman entropy. 



5.4. Topological phase in strategy space for 
multi-player game 

We study the topological quantity in the strategy space 
of n-player game in this section. The total payoff function 
of the n players could be view ed as Hamiltonian matrix. 

H = {mi,u 2 , ...,«,■}. (5.15) 

A strategy vector of the many players's strategy space 
reads 

. = K«.f...ft (5.16) 

s£ n ^ is the A:th strategy of the nth player. u n maps this 
strategy vector into the payoff value of the nth player, 
i.e., u n s = u n . When the Hamiltonian matrix of the 
payoff function operates on this vector, it produces the 
eigenvalues 

Hs = {lii,1i 2 ,...,ttn}- (5-17) 

The players fight against each other to reduce the differ- 
ence between their payoffs 

Ai = ui - u 2 , A 2 = u 2 - u 3 , A n = u n - ui. (5.18) 

In a quantum system, u p could be interpreted as the pth 
energy level. The Nash equilibrium solution is governed 
by the coexistence equation 

D(A/s) = {Ai,A 2 ,...,A n } = 0. (5.19) 

The war between energy levels breaks out in the degen- 
erated eigenspace. The players take their strategy to de- 
crease the energy gap between them. Those points at 
which the energy gap vanishes are the core center of vor- 
tex, they are the energy level crossing point. The ith 
player choose strategies according to other players' strat- 
egy, thus his strategy vector is a map from the other play- 
ers's eigen-strategy to his own space ipiil)- I n analogy 
with conventional Berry phase of quantum mechanics, we 
proposed a general topological quantity (See Appendix[j]). 
The topological phase of many players reads 

fi< = ie t j.. .kiidjMl^djjMl))--- (d lk Ml)\dj, Ml))- 

fli is actually the Ricmannian curvature tensor in the 
strategy vector bundle. This quantity is equivalent to a 
topological current of the payoff functions 

fif = JJ(5(Aj - Aj)<2Ai A dA 2 A ■ • -dA„_i A dA n . 

From our general generalization of conventional Berry 
phase, the most general topological quantity is the Chern 
character. In this strategy vector bundle, the Chern char- 
acter can be derived from the exponential map, 

Ch{x) = Exp(d{il>(~f)\ A d|V>(7))). (5.20) 
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These topological quantities are good candidate to mea- 
sure the entanglement between different players in strat- 
egy space. One can see from Eq. (|5.20j) . topological 
charges focused on the equilibrium solutions where all 
players get the same payoff, there is no difference among 
any two players. But remember these equilibrium point 
are isolated point, they are the center of vortex. A tiny 
deviation from these point would break the equilibrium. 

6. GAME THEORY OF MANY BODY PHYSICS 

6.1. Cooperative game and classical many body 
system 

The particles of a physical is the players of a game. 
They act following the least principal, the aim of the 
players is to decrease the total energy. Some of them 
may united to form a subgroup due to local potential. 
Two subgroups may fuse into one when this fusion can 
make both of them more stable. It will be shown the 
many player cooperative game(see Appendix |F|) is just a 
physics system. 

We consider a 3-person cooperative game {1,2,3}. 
We choose the payoff function of this game as 
the conventional Hamiltonian, the first step of the 
cooperative game is to find out all subgroups 
{0, 1, 2, 3, (1, 2), (1, 3), (3, 2), (1, 2, 3)}. The energy fane- 
tion maps each subgroup to a real number, 

£0 = 0, El, £2, £3, £\2, £23, ^13, £l23- (6-1) 

This energy mapping may have various formulas. The 
total energy of the three particle is E{N) = s\ + e% + 
£3 + £12 + £13 + £23 + £123 ■ The super-additivity requires 
£2 + £3 < £23- e(23) = £23 — £2 — £3 is the difference of 
payoffs called excess. In fact, this is a nuclear reaction in 
physics. In game theory, people are looking for optimal 
self-content energy vector of the n particles e* and the 
excess value pair . In the eyes of physicist, this is group 
expansion algorithm. 

For a standard Hamiltonian system, 

- 1 

H = J2^P 2 i+V(qi,...,q N ), (6.2) 

i=l 

we take N particles as players, their position qi and mo- 
mentum P l as their strategy, the energy is the payoff. 
The target of this game is to find the stable ground states. 

To investigate the phase transition of this game in 
the frame work of our topological phase transition the- 
ory, one only need to replace the output function O u t (7) 
with the Hamiltonian H, or any other independent in- 
tegrals in the complete set (<f>i — H, 02, (f> n ) which 
are in involution. The strategy vectors is taken as 
7 = ((ft, QjViPii ■■■,Pn)- Then all the conclusion on 
topological phase transition we obtained in the previous 
section holds here. 



However the physical system we encountered always 
have a very larger number of particles, usually N ~ 10 23 . 
It is impractical to control the position and momentum 
of particles. Therefore we take a different way to model 
the many body system. We focus on a few interaction pa- 
rameters and integrate the configuration space out. For 
example, we choose the Hamiltonian, 

N 1 

H = J2 2 P * + ^ V i(<H, -' qN ^ + T^fai, Qn) + 

i=l 

The interaction potential relies on n players above con- 
figuration space, 

y(g)= 7 V 1+7 V 2 + ... + 7 »y„. (6.3) 

where {71, 72,...,} are physical parameters, and 
Vi = Vl(qi, gjy) is the potential functions. We 
integrate out all the uncontrollable position and mo- 
mentum variables in the Helmoltz free energy on 
configuration space, then the effective output function 
for a finite player game is derived, ^(71,72,...) = 
-(2(3)- 1 Log(n/(3) - /(/?, 7l , l2 , ...)//?. ~ with 

f (0,71,72, ■■■) = j? L °9 J ' d N q exp[-0V(^,'Yi,j2, ■■■)}■ 

The free energy function is merely one special compo- 
nent of output vector field corresponding to Hamiltonian 
function. We can choose any output function O u t on 
the strategy manifold, and choose a arbitrary component 
of the compete set whose members are in involution as 
Hamiltonian function <pu, then we introduce the renor- 
malization group transformation upon the output vector 
field. The basic tangent vector field cf> to characterize the 
topology of the extended strategy space is given by the 
variational principal 

S p [UOut(7i,72,-)} = J d N qd N p exp[-/?0 g ]ffl[?7O Mf ] 
S0P N Jd N qd N p exp[-f3<j> H ] ' 

A special subset of renormalization group transformation 
is Lie group U = e l6L . Now we can apply the topologi- 
cal phase transition theory to study the phase transition 
on the strategy space extended by the interaction pa- 
rameters 7 = (71,72,...). The conventional conjugate 
momentum and position variables have been integrated 
out, they are out of the set of players for the phase tran- 
sition game. Of course, people can bring them in the 
game, in that case, the number of players is too large to 
manipulate, it does not tell us any practical information. 

Generally a general dynamic system can be viewed as 
a non-cooperative game between players 71, 72,... "f n - 
If we view the classical system as a coopera- 
tive game, we can divided the effective interac- 
tion parameters {71, 72, 7n } into several groups 
{(71,72), (•••7i), — >7n}i an d consider the interaction be- 
tween these groups. Then we analyze their coalition, 
their excess value pair, their payoff, etc. We could get 
the most effective information about the system through 
the phase coexistence equation. 
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6.2. Topological phase transition of quantum many 
body system 
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The topological phase transition theory in this paper 
aims at the most general systems, no matter they are 
classical physical systems or quantum physical systems, 
or biological system, or social system. For a quantum 
many body system, any output that responses corre- 
sponding to certain input can be used to detect a phase 
transition. The output O u t could be any physical observ- 
able, such as the correction to ground state energy SE g , 
the correction to thermal potential <5f2, the single-particle 
current operator J = J d 3 xj, the number density oper- 
ator < p(x) >, the spin density operator < a{x) >, the 
free energy F, susceptibility x, specific heat Cjj, corre- 
lation length f , compressibility kt, ■ • • , and so on. The 
input are any physical parameters, such as magnetic field, 
electric field, radiation, temperature, pressure, on-site re- 
pulsive interaction, chemical potential, density, potential, 
scattering length, neutral current, electron wave, probing 
laser beams, enzymes, and so forth. Our topological 
phase transition theory can be successfully applied to ex- 
plain the quantum tunnelling of the magnetization vector 
in ferromagnetic nano-particles 22j . 

We take a specific model to demonstrate the appli- 
cation in the following. A gas of ultracold atoms in 
an optical lattice has provided us a very good exper- 
imental observation of superfluid-Mott-insulator phase 
transition The system is described by the Bose- 
Hubbard model 



H 



E 



(6.4) 



where the sum in the first term of the right-hand side 
is restricted to nearest neighbors and a\ and dj are the 
creation and annihilation operators of an atom at site i 
respectively. J is the hopping parameter. U corresponds 
to the on site repulsion between atoms, and \x is the chem- 
ical potential. This Hamiltonian admits two conflicting 
forces, J and U, player J drive the particles hoping from 
one site to another, but U repulse any particles jumps 
to his sits, this prevent the particles from moving site 
by site. The quantum phase transition is governed by 
the two players, when the hoping player is dominated, 
the ultralcold atoms can easily hop in the optical lattice, 
this is the superfluid phase. When the on-site repulsive 
force is dominated, the particles is strongly repulsed by 
its neighbors, so it would be difficult for him to move to 
other house, then the ultralcold atoms have to stay at 
home, this is the Mott-insulator phase. 

Using mean-field approaches, the ground state energy 
of Bose-Hubbard model can be generally expressed as 
the functional of J and U, i.e., E g = E g (si,J,U). 
Following perturbation theory up to the second order [23j. 
the variation of ground state energy 5E g is 
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FIG. 6: The phase diagram of the Bose-Hubbard Hamiltonian 
for g=l. Inside the curve is the Mott-insulating phase, outside 
the curve is the superfluid phase. 



From our definition of the first order phase transition, the 
boundary between the superfluid and the Mott insulator 
phases should be decided by the equation 5E g = 0. We 
have plotted the phase diagram(FIG. 6). This phase di- 
agram is in agrement with the familiar phase diagram of 
superfluid-Mott-insulator phase transition; 23j. It is also 
the same as the phase boundary obtained from the min- 
imizing the free energy 2J| , where they take the hopping 
term as perturbation. 

In this phase diagram, each phase is assigned with a 
winding number. The winding number /3kVk — Wfc is 
the generalization of the Morse index in Morse theory. 
Its absolute value /3k measures the strength of the phase. 
t]k =sign{</> , <fi 2 } = ±1, represents different phases. In 
the dark area circulated by the curve, r)k =sign{0 , <p 2 } = 
+1, it represents the Mott-Insulator phase, while outside 
the curve, r\k =sign{<^ 1 , (f) 2 } = —1, it represents the su- 
perfluid phase. On the curve, {tp 1 , </> 2 } = 0, that is where 
a quantum phase transition takes place. 

Most of the recent phase diagram of the quantum phase 
transition are focused on the first order. As for the pth- 
order quantum phase transition, it is characterized by 
a discontinuity in the pth derivative of the variation of 
ground state energy. The order parameter field of the 
pth order Quantum phase transition can be chosen as 
the (p-l)th derivative of 5E g with respect to U and J, 
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dfj 1 SE g . The equation of 



coexistence curve is given by {(j) 1 , </) 2 } = D(-) = 0, 
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here we have taken 71 = U, 72 = J. From this equation, 
one can arrive the phase diagram of the pth order quan- 
tum phase transition from the equation above, and find 
some new quantum phases. 

The above discussion is based on the correction to 
ground state energy, now let's choose the output field 
as thermal potential, O u t — A51. One can do some per- 
turbation calculation up to the second order, and derive 
the Green function, 



G-\p, l c n ) = e - e 2 S> + 1){ - na -gg±* } (6.7) 
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The two parameters are {/i.U}. The correction of the 
thermodynamic potential is 

Afi = fl — fi = -s^H-pG^faioJn)]. (6.8) 

The first order phase transition is given by Af2 = 0. 
Put Eq. dlTT]) into Aft = 0, we have G _1 = -T. At 
zero temperature, T = 0, it leads to G _1 = 0. This 
is exactly the condition for the superfluid-mott-insulator 
phase transition in Ref. t 25j. 

Now we see different physical observables are only dif- 
ferent outputs, they lead to the same phase boundary. 

6.3. Elementary excitations and momentum space 

One of the best merits of quantum many body theory 
on lattice is we do not need to consider canonical po- 
sition coordination and its conjugate canonical momen- 
tum. The wave vector has only three spatial components. 
Although we can not manipulate them, they are very 
good candidates for the player of a game. 

The inverse of the wave vector is wave length which 
is quantized due to boundaries in a finite solid state lat- 
tice. Like the string fixed at both ends, the wave vector 
of the standing wave modes only take a few discrete val- 
ues. The standing wave patterns of the electron wave 
and atomic wave in lattice is viewed as elementary exci- 
tations, or quasi-particles. If we think the electrons as 
soldiers, the elementary excitation is the collective danc- 
ing of millions of soldiers. When the soldiers are confined 
in a square battlefield, they self-reorganized into certain 
dancing modes to avoid crashing each other. In fact, the 
independent spatial components of wave vector are their 
commanders, this is a war game between wave vectors. 
These oscillation modes, or elementary excitations, are 
the surviving strategies of the commanders in this bat- 
tle. 

In quantum many body system, the elementary exci- 
tations sit at the singular points of the green function. 
The external response can be derived from green func- 
tion by linear response theory (see Appendix |K|) . These 
responses are the output vector field of a quantum many 
body game. For the most general case, we set the 
Green function on a m dimensional momentum space, 
P = (Po = E,pi,p2,...,p m ). In analogy the topologi- 
cal current of Riemannian curvature tensor on m dimen- 
sional manifold, we introduce the topological current of 
TV-point Green function G in momentum space, 

B G = e tJ ... kl (d Pt G\d P] G)...(d Pk G\d Pl G), (6.9) 

where we denote \d Pl G) — d Pl G(p) and (d Pk G\ — 
d Pl G\p). The inner product between Bra (|and Ket |) is 
an integral, i.e., (d Pt G\d Pj G) = J dpd Pt G^ N (p)d Pj Gi(p). 
Use this topological current, we can construct a more 
general topological action 



This topological action measures the non-trivial topology 
of Green function which maps the momentum space to 
an external response. The Green function correspond- 
ing to the exactly solvable terms of hamiltonian can 
be exactly derived, the non-exactly solvable part is al- 
ways calculated using perturbation theory. According 
to Dyson equation, the difference between the inverse 
of exact Green function and the exactly solvable Green 
function gave us the correction to energy, 

SE = G{p, E)- 1 - G Q (P, Ey 1 = E(p, E). (6.11) 

usually people call it self-energy T,(p,E). In this case, 
the nontrivial topological current of the N-point Green 
function focus on 

B p = Tr[d{S\ A d|E) • • • Ad(E| A d|E) 2p ], (6.12) 

again we have the holographic topological action 

Ch(M)v = Tr(exp -d(£(p, E)\ A d|E(p, E))). (6.13) 

The self-energy here is the output vector field. The mo- 
mentum vector and frequency are the strategies. The 
surviving strategies are where the stable elementary exci- 
tations arise. One may apply a Lie group transformation 
to self-energy to check if there the discontinuity exist. 
Usually we take 50(4) whose element can be written as 

U(9) = e eL = Y2 7j( Ld Y\ where L is the generator of 
50(4), it is operator constructed from the momentum 
vector P and energy E. If there is discontinuity, we in- 
troduce the vector field <f — 5u ( e ^P> E } m Then we can 
use the phase coexistence equation {<//, cjp , <f> k } = to 
decide the phase boundary. 

We always assume the periodic boundary condition in 
solid state physics, this boundary condition created a 
torus lattice manifold. This naturally leads to a topo- 
logical constrain in momentum space. If we take the 
components of momentum as the player of a game, their 
strategy space is the momentum space. If the momentum 
space is noncompact, there is no topological constrain on 
strategies. If the momentum space is a compact mani- 
fold like sphere or torus, there is a nontrivial topological 
constrain on strategy manifold. 

To avoid the non-figurative reasoning, we consider the 
superconducting pairing Hamiltonian [261 ] . 

He = ^(^-4^ + AyC^^ + A* jCji c lt ) - fl^ cJ ff c i(J . 

ij i(T 

(6.14) 

If the system is translationally invariant, then tfj = t a (i— 
j) and Ay = A(i — j). We further require that hopping 
rate for the up-spin and down spin is the same, t\- = t^, 
then the up-spin and down-spin have the same kinetic 
energy (k) = ^(k) = £(fc) = £• e^^ij) - fi. The 
momentum representation of Hamiltonian H e is simpli- 
fied as 



Ch(M) = Tr(exp -d(G(p)\ A d\G(p))). (6.10) 



H e (k) = R(k) ■ a, (6.15) 
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where a = (<J X , a y , a z ) are the Pauli matrices, and R = 
(R x ,R y ,R z ) = [ReA(k),-ImA(k),Z(k)] with A(fc) = 
J2j e~ tk ' T3 ' A(j) as the gap function. Using the Bo- 

goliubov transformation, ctk = UkCk — VkC_,, ar k = 
u k c l ~ v k c -k and considering [ak,H e ] — E^a^ for all k, 
the effective Hamiltonian becomes H e = ^ fc EkCt\.ctk + 
const with Ef. > 0. The effective energy of quasi- 
particle is Ek = + |Afc| 2 . The Anderson's pseu- 

dospin vector [27| 

_ (ReA k , -JmA fc , 6) 

u= = . (6.16) 

The pseudospin vector n describes a mapping from k- 
space to a sphere S 2 in the pseudospin space [28j. This 
map defined a topological number to classify the topology 
of the momentum space. 

As observed in the calculated results for the berry 
phase curvature in SrRu03. ft has a very sharp peak 
near the T-point and ridges along the diagonals [29] . The 
origin for this sharp structure is the degeneracy and band 
crossing. The transverse conductivity a xy is given bv[30| 

k " r 

We applied the topological current theory [l2| to the mo- 
mentum space of this pairing Hamiltonian, and proved 
that the curvature only exist at the solutions of R = 0, 
i.e., 

n xy = J2 w iW-K) ( 6 - 18 ) 

k 

While W\ is just the winding number of R around the 
Z-th solution R(ki) = 0. 

If we take the whole two dimensional system as a 
war, the soldiers are the two momentum components 
k = (pi,P2)- Each pi has a strategy space which is 
his available value across the whole momentum space. 
Pi fights against each other to steer the navigation of 
quasi-particle. If the whole system is homogeneous and 
external potential is not anisotropic, pi and P2 are ex- 
changeable. In real space, the electrons in a homogeneous 
system has no bias direction, they circular around each 
other. Two dimensional system is rather special, the elec- 
trons have to dance around each other to avoid crashing. 
When the external magnetic field is present, the elec- 
trons change their pattern of motion in real space, they 
are adjusting their wavelength and oscillation modes in 
the mean time. The electrons are spontaneous organized 
into collective oscillating modes in which any two of them 
form pairs. 

The interference of the electron wave forms periodic 
distribution in the two dimensional momentum plane. 
Generally speaking, p\ and pi are much like coworkers, 
they balance their inner conflicting profits to fit the envi- 
ronment. The best surviving strategy focused on the ex- 
tremal point of the external potential, at these points, the 



quasi-particle or collective oscillation modes do not cost 
energy any more. Those solutions are given by R = 0. 
These periodically distributed extremal points are the 
basic channels through which the electron like to travel. 
Each of these channels may be view as the center of cy- 
clone, they carry a winding number. The sum of the 
winding numbers of all these channels is under a topo- 
logical constrain which comes form the topology of mo- 
mentum manifold. 

The topology of momentum manifold strongly depend 
on boundary condition. It seems strange for a brick of 
material with a rough surface to have some periodic or 
other boundary condition for the wave function of elec- 
trons. In fact, like the droplet of water in vacuum without 
gravitational field, their surface are spontaneously com- 
pact sphere. Compact boundary is one way for a system 
to fit the environment under certain environment. In 
other worlds, compact boundary condition is the result 
of evolution following the Principal of Least Action(The 
principal we talked bout here is the generalized principal 
of least action in the first several sections of this paper). 

Therefore he electrons are organized into certain com- 
pact boundary condition by themselves in order to save 
energy or to form a more stable state. The topology of 
their momentum space is modified according to exter- 
nal inputs, such as chemical potential, magnetic field, et 
al. Anyway, if the Hamiltonian relies on many physi- 
cal parameter — 1,2,...}, i.e., H = H (7 1 , 7 2 , ...), 
the coexisting boundary at which they reach a balance is 
given by 

D(i>/^) = {$i,4> 2 ,...,$ n ]=0. (6.19) 

If one of the operator is taken as Hamiltonian itself, <j> = 
H, this is nothing but the complete set of dynamic system 
in parameter space. 

6.4. Quantum many body theory of game 

The player of war consists of a large number of soldiers 
which act like identical particles. The soldiers split into 
several subgroups, inside every subgroup there is a com- 
mon agreement constituted and obeyed by all the mem- 
bers of the association. Thus we may take mean field 
approach inside the subgroup whose member move in a 
mean-field potential formed by all the other members. 

Generally speaking, the subgroup consists of angel 
players is stronger than that of devil players. Angle play- 
ers are bosonic players. The collective strategy is the 
symmetric combination of the personal strategy, 

|5 Bos (l,2,...,n)) = yE[F(| S W)| S ( 2 ))....| S («))). 

The devil payer are fermionic players. Their collective 
strategy is a totally anti-symmetric strategy, which can 
be written as Slater determinant, 

\ S ^r) = E ^^...^ Qia2 - Q "(l4 ) )l^)--I^ ) ))- 
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The above discussion is only suitable for some oversim- 
plified system. When it comes to a general multi player 
game, we do not distinguish which one is angel player and 
which one is devil player. They are all rational players. In 
every round of the game, the player choose one strategy 
from his strategy space, so do other players. When they 
check their payoffs, what someone win is what someone 
lose. We may call the individual loser an angel player, 
and call the individual winner a devil player. 

In fact, the players share the same strategy space. A 
wining player takes good strategies, a losing player takes 
those bad ones. It does not make any sense to talk about 
the good or the bad of the strategy itself, the same strat- 
egy may bring different results in different situations for 
a certain layer. Every combinatorial series of strategies 
corresponds to a payoff distribution to individual play- 
ers. If we take particles as player, a physics system always 
contain millions of particles, it is impossible to calculate 
every particle's payoff. That is why we need statistics 
physics through which we summarize the information of 
payoffs into a few statistical observables, but we lost the 
information on individual payoffs. We choose a more 
practical way to grab the key information of many parti- 
cle system. The player of game is taken as the different 
interactions that govern the microscopic behavior of par- 
ticle. 

A vacuum state |0) is defined by a state that no player 
confirm his strategy. When the ith player has chosen 
his strategy \sal), and place his card on the table, it 
means that a strategy s a is generated from vacuum, i.e., 
|sq' ) ) = sj Q |0). One player's strategy is not isolated, it 
is always entangled with other players 's strategy, since 
whenever he makes the decision he must reflects other 
player's choice, the collective strategy is highly entangled 
states. The anti-strategy of the ith player comes from the 
combination strategy he received from the rest players, 
we denote it as 

(S % \ = (0|F(s Ql ,s Q2 , ...,s ak , ...s a „), (k^i). 

When the ith player encounter the jth player, his 
payoff is Uij — Tr^j [(5*|sq )]. We can further define 
the density matrix p = YljPj\Sj)(Sj\> an( A introduce 
the familiar Von Neumann entropy to measure the 
entanglement. 

Pairing mechanism in the war game 

As shown in the prisoner dilemma, two altruistic 
players can get their best payoff without communication 
and cooperate, they trust each other and both choose 
the strategy for the welfare of other people. But two 
selfish player can not get the best payoffs unless they 
communicate and cooperate, in order to prevent the 
other player's betray, they would like to bound into 
a pair so that any betray would damage himself, this 
increases the entanglement between them, in physics, 
two spin in magnetic field is a perfect physical system to 
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FIG. 7: Two spins in magnetic field is a good demonstration 
of prisoner dilemma in physics. 

realize a prisoner dilemma. A direct physical observation 
is that the particles with spin is lean to form pairs, and 
the magnetic field would strengthen the entanglement 
between two particles with opposite spins. 

A physical system to demonstrate the Prisoner's 
dilemma is two spins in magnetic field. Two players are 
the two particles: Alice and Bob, they both carry mag- 
netic dipole momentum ji. The strategies: confession and 
accusation corresponds to spin up and spin down, the ex- 
ternal magnetic field is external law. Their loss function 
is their energy in this system. They fight to minimize 
their loss function. The payoff matrix is 
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A is the energy shift due to the interaction between the 
two particles. This is the classical picture of two-particle 
game. There was experimental realization of prisoner 
dilemma on nuclear magnetic resonance quantum com- 
puter (3lj |. 

A general n-player prisoner dilemma may be summa- 
rized into a n-spin Hamiltonian. A player has two strate- 
gies: spin-up and spin-down. The payoff function is the 
Hamiltonian, 

H = J2hi<ri + J ij Y,44- (6-20) 

i ij 

If we take the local coupling J y = J and consider only 
the nearest neighbor interaction, this Hamiltonian be- 
comes traditional Ising model. The payoff matrix of this 
n-player prisoner dilemma game for Ising model is diag- 
onal block. A group of soldiers in a war also encounter 
prisoner dilemma. Two soldiers helping each other is 
stronger than the two that they do not help each other. 

We consider a war on two dimensional lattice. The sol- 
diers are particles, they are divided into two equal groups, 
A and B. The particles of A is the anti-particle of B and 
vice versa. In the beginning, A and B are separated two 
parts on the lattice. When the war break out, they rush 
to each other and began the combat. A particle only 
fight against his nearest antiparticles and help his near- 
est friend neighbors. If the soldiers are identical particles 
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at the same level, the interaction between particles does 
not go too far, it is confined in a local region. If we 
assume all the soldiers only fight against enemy soldiers 
for his own survival, they do help friend neighbors, the 
dynamic process only includes 



{£ij §\j Bj , £jiS J B .s Ai , AieA,BjeB}, (6.21) 

where s represents the soldiers. This is a system com- 
posed of free particle and antiparticlcs. In fact, some 
soldiers of G may help each other, they grouped into 
pairs to fight. So does the soldiers of G. Then we have 
to consider the pair interactions. The war now contains 
three solders interaction and four soldiers interaction, 

£ ij 4* B i' £ji S Bj S Ai , 



V 3 s A / A J B] , ... 
V 4 s A s A s B J Bl 



, Sji S B . S Bk S Ai , 

y4 s Bi s B .s Ak s Al 



(6.22) 



If we recall the BCS pairing mechanism in superconduc- 
tor theory, the war between electrons does not have three- 
particle interaction. Even two electrons separated far 
from each other, they can act as one pair. They do not 
even know each other, since identical particle arc indistin- 
guishable. This interesting phenomena in the quantum 
world may found some naive analogy in war. During the 
war game, we can assume a pair of soldier are confined 
in the lattice, they kept running from one place to an- 
ther, but helping each other all the time. They shot any 
enemy that try to kill his partner. This is a long range 
interaction. 

There are two types of pairs, angle pair and devil pair. 
The angle pair is more powerful than the devil pair as 
a whole. They love each other and trust each other, so 
they can separated from each other in a longer distance. 
But devil player have to stay closer to prevent his partner 
from betraying him for his own survival, their correlation 
length is shorter. 

There are hierarchical structure in army. The soldiers 
are organized into different groups in which a large num- 
ber of them fight as a whole. These hierarchical struc- 
ture are well kept in the region far from the coexistence 
boundary, but when combat begins in one battlefield, 
all these hierarchical structure becomes meaningless, the 
general or the captain at high plays no different role as 
an ordinary soldier at the lowest level, they are just the 
same individual fighters with arm. It is the hierarchi- 
cal structure or inner structure of a system that differs 
it from others. When the phase transition occurs, all 
these hierarchical structure are broken, men are born to 
equal, men are also died to equal. This is the origin of 
universality class. 

For a further consideration of war game, we have to 
take into count of different hierarchical structures. The 
commanders are the critical players, they are dressed up 
by a group soldiers, and fight as a whole around the bat- 



tlefield. We may take Kadanoff block procedure to sum- 
marize all the many player interaction to the commander. 
The many body interactions, such as 



v n 4 ... a 



SA k 



SA n 



are renormalized into partition function. The war is go- 
ing on between the renormalized commanders, S B and 
S A , where 



(4 -4) 



S B , (4 ...s A )^S A . (6.23) 



Again we may consider the pairing interactions among 
these commanders, 



i £ ji "4 4t ^Ai , 



v s 3 s A s Ak s Bj , 
v s s A . s A _ s Bk s Bl , v 4 4, 4 4 4 • 



(6.24) 

Repeating this renormalization group transformation, we 
transform the sophisticate many body problems into 
something we can handle. This is the fundamental spirit 
of various numerical renormalization group method de- 
veloped in quantum many body system. 

The phase boundary appears when the fighting groups 
reach a balance at certain frontiers regions. Friend and 
enemy are maximally entangled in the coexistence state. 
One can not tell who is who. In order to study the sym- 
metry transformation at the coexistence boundary, we 
introduce the position displacement operator and its con- 
jugate momentum operator, 



Jt 



V2 



V2 



(6.25) 



Then the angular momentum operator is given by Ljj = 
(XiPj — XjPi). The generator of the symmetry trans- 
formation may be denoted as L — i ^ . OLij , through 
which we derived the Lie group element 



U 



„iOL 



-iJLiiOiSlst-SiSj) 



(6.26) 



The quantum output vector O ut (Si, Sj, ...) is expressed 
by the operator of renormalized soldiers. Then we can 
detect the phase transition by the transformation 



U( P )OutU^ = e l9L 6 ut e- 



ie-L 



(6.27) 



This transformation equation determined the evolution 
of the quantum output in the vicinity of phase coexis- 
tence region. Inside one phase, or inside one army, this 
transformation is continuous up to infinite order. When 
never the symmetry lost, an enemy popped out and ini- 
tiated the revolution. 

The uncertainty principal shows up across the whole 
phase diagram. In the phase diagram, the war game is in 
different squeeze state for different stage of combat. The 
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transformation operator (|6.26l) is a special squeeze oper- 
ator. A general squeeze state[37J generated from vacuum 
is 

|V») = e- i S«W-«*^i)|0). (6.28) 

When the two players began rushing toward each other 
but still do not touch each other, their fluctuation of 
spacial position is much higher than that of their force. 
When they meet in the phase coexistence boundary, the 
fluctuation of spatial position is around the boundary, it 
is much lower than the conflicting force. At the Nash 
equilibrium state, the phase boundary is almost sta- 
tionary, but both of the two players summon up their 
maximal force to against each other, they reach a dan- 
gerous equilibrium state on the boundary. This is an 
extremal squeezed state, the fluctuation of position is 
highly squeezed, but the fluctuation of conflicting force 
is very strong. 

So we may take two Hermitian operators 7 and P 7 , 7 
denotes the position in phase diagram, P 7 represents the 
corresponding force. The commutation relation derive 
from the second quantization is [7, P 7 ] = iC. We can 
calculate the uncertainty of an operator A by (AA) 2 = 
(ip\A 2 \tp) — ((ip\A\ip)) 2 . According to Heisenberg uncer- 
tainty principal, A7AP 7 > h\C\. The uncertainty of the 
position operator is 

A 7 «i|C| (6.29) 

at the coexistence state. While the uncertainty of corre- 
sponding momentum is much stronger AP 7 3> h\C\. 

The angel pairing and devil pairing are two stable 
squeezed states. The angle pairing has less uncertainty 
about their trustworthiness, so they have larger uncer- 
tainty on their spatial connection. But the devil pair 
has larger uncertainty on his partner's trustworthiness, 
so they strengthen their spatial connection, this reduced 
their spatial uncertainty. No matter which kind of pair- 
ing it is, the total information in the two pairing states 
should be conserved. If we can make sure a system doe 
not lose information, we can maintain the pairing to a 
stable level. Superconductor is such a system, below the 
critical temperature, no chaos invades into it, no infor- 
mation escapes out of it. There is a constant electric 
flow. 

According to quantum field theory, symmetry means 
conserved quantity, a loss of symmetry indicates a loss 
of conserved quantity. What the old phase lost becomes 
the generator of the new phase. The squeeze transforma- 
tion on output field is a way to detect the loss of infor- 
mation. Imbalance is the original engine of development. 
For physicist, the physical measurement always affect the 
output state itself. We can not detect a state without 
changing it, this minor change is neglectable when it is 
done in a stable phase. But at critical region, the mi- 
nor effect of detection might induce significant effect. In 
other words, when a war is going on, once the spy step 
into the frontier, he is force to fight for his own survival, 
it would be too hard for him to sent out any information. 



7. SUMMARY 

(l)game theory of general phase transition and 
Renormaization group transformation 

Phase transition is a much more universal phe- 
nomena than the conventional phase transitions in 
physics. Similar sudden changes arise from evolution 
of biomolecule, self-organization in nanoscale systems, 
cosmic evolution, and so forth. The general phase 
transition breaking the envelop of physics can be viewed 
as a non-cooperative game of many players. The players 
try different strategies to win. Whenever a winning 
player fails and becomes a loser, a phase transition 
would occur. 

When the players represent different interactions be- 
tween the elements of a complex system, a winning player 
is a dominant interaction which governs one stable phase 
of the system, the other phases represented by loser's in- 
teraction are suppressed. When we tune the interaction 
parameters, we are changing the strategies of the players, 
the winning interaction may becomes weaker and weaker. 
There is a critical point at which the losing players grows 
stronger enough to balance the previous winning player, 
then we arrived the Nash equilibrium point. At this point 
, no one wins, and no one losses, it is the coexistence 
point of the new phase and old phase. The Nash equilib- 
rium point is an unstable saddle point, any tiny deviation 
would decided the fate of the coexisting new phase and 
old phase, there is only one winner left to dominate the 
phase structure of the system. 

Usually the Nash equilibrium of a war game is not de- 
rived in one round of combat. The players have to play 
many rounds of negotiation to find the optimal strategy. 
The renormalization group transformation theory is ac- 
tually one theoretical description of this kind of game 
process in physics. In physics, the players are the inter- 
action parameter or physical parameter. For example, 
we take the Ising model as a game, the two players are 
the spin coupling interaction and external magnetic field. 
The stable phase is determined by the dominant interac- 
tion. 

At the first round of renormalization group transfor- 
mation, the two player accomplished the first round of 
combat, they made an agreement on the most part of 
boundary between their domains. But the unsettled part 
is occupied by the winner. So the loser initiated the sec- 
ond round of renormalization group transformation to 
get some of his losing back. As this transformation goes 
on, the unsettled boundary between the players becomes 
less and less. Finally, they reach the Nash equilibrium 
point, that is the fixed point of the renormalization group 
transformation. This point is where the phase transition 
occurs. 

Usually to obtain the exact critical point, one has to 
perform the renormalization group transformation to 
infinite order. However, it is too far to reach in reality, 
we always cut it off at certain order. It inevitably gives 
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birth to a loser and a winner to certain order. Although 
the players at that order maybe are fighting for only 
one thousandth of a penny, but the loser is loser, winner 
is winner, the title is fixed. The order of cutoff is the 
order of phase transition, we can take renormalization 
group transformation on output function to observe 
different order of phase transition. When we apply 
this definition to thermodynamics physics, it naturally 
unified Ehrcnfcst definition. 

(2)Universal coexistence equation and 
topological conjecture of scaling laws in 
fractal space 

The game is played on the strategy base manifold, 
so a phase transition is related to the topology of the 
strategy base manifold. If the strategy space is noncom- 
pact, we do not need to consider topology at all, in that 
cases, the phase transition governed by game dynamics. 
We may call it a non-topological phase transition. When 
the strategy manifold is compact, there would be a 
topological constrain on the player's strategies, then the 
topological effect would appear in phase transition, we 
call it a topological phase transition. 

The topology is intimately related to the fixed points 
of renormalization group transformation flow on strat- 
egy manifold. We introduced the general flow vector 
field, and showed that the sum of the winding number 
around the surviving strategies is a topological number. 
These surviving strategies with opposite charges annihi- 
late at the universal coexistence curve, whose equation is 
{(f) 1 , (f> 2 , (j) 2 } = 0, <j) 1 are the component of the flow vec- 
tor field of output field. In the thermodynamics physics, 
the output field may be chosen as free energy, the players 
are temperature and pressure, one can verify that this co- 
existence equation unified all the coexistence equations 
at different order classical phase transition. The phases 
separated by the coexistence curve are assigned with dif- 
ferent winding number of opposite sign. 

The strategy space is a fractal dimensional space. It 
is easy to see this point from a war game. We focus 
on the battlefield, it is a war between two armies. If 
we look closer, it is the war between corps. We may 
continue the magnify the battle field, one would see, it is 
a war at all different scales, it is the combat between two 
individual soldiers at any local region. A war game may 
be summarized into the war between two commanders 
by Kadanoff-block procedure. There is an intrinsic 
fractal structure around the critical point. So we make 
the hypothesis that a local neighborhood on the strategy 
base manifold is homeomorphic to a fractal dimensional 
space. Then the output field may be expanded by the 
polynomial of power functions with fractal dimension 
index. In analogy with the observables defined in 
the statistical mechanics, we introduce the tangent 
vector field of the output field as observables. These 
tangent vector field may be approximated by fractal 
polynomials in the vicinity of Nash equilibrium solution. 



By substituting these fractal polynomials into the uni- 
versal coexistence equation, one can derive the scaling 
relations. As we shown, the coexistence operator has a 
degenerated subspace, the coexistence equation bear a 
topological origin, so these scaling relations are universal. 

(3)Many body physics of games 

Game theory is actually one different way to see 
many body system. For example, a collection of millions 
of ultralcold atoms trapped in optical lattice can be 
described by Bose-Hubbard model. There are two key 
parameter, the hopping parameter J and the on site 
repulsive interaction U. We can model this system 
as a game between U and J. Player J command the 
particles hoping from one site to another, drive them 
into superfluid phase. While U direct the soldiers to 
stick to the lattice site. If J wins, the ultracold atoms 
becomes superfluid. If U wins, the ultracold atoms arc 
confined in the lattice, they form Mott-insulator phase. 
The transition point is the Nash equilibrium point of 
this game. The two phases coexist at the critical point. 
Another way to model the millions of ultracold atoms is 
to take every atom as a player, and their states at lattice 
sites are strategy space, this reproduced the standard 
quantum many body theory. 

A practical approach to a given many body system is 
first to find out the main different interactions that are 
players of the game. Secondly, we choose some output 
quantity to measure the states. Physicist usually take 
the external response function, these output field are the 
payoff function of the game. The tangent vector field es- 
tablished on the hypersurface of these response function 
is the fundamental vector field to study the topological 
effect of the strategy space. Applying the universal co- 
existence equation, we can find out the basic structure 
of different phases. This scheme holds both for classical 
systems and quantum systems. 

Game theory only provide us a mathematical frame 
work to understanding the general property of a system. 
The topological phase transition theory developed in this 
paper is a general mathematical results, it does not de- 
pend on any specific model in physics, biology or social 
system. When it comes to some specific physical system, 
we may apply the conventional theory for that special 
system to derive the specific relation between input pa- 
rameters and output field, or one may directly conduct 
some experiments to find some approximated relation. 
As long as one derived the specific relation between out- 
put and input, just substitute it into the universal coex- 
istence equation, one can get basic phase diagram. 

In fact, quantum many body theory may provide us 
new understanding to many player games. We developed 
the density matrix theory of many player game, it was 
shown there is a one-to-one correspondence between fixed 
point of many player game and the self-character vector 
of the player. Every quantum phase may be represented 
by a state function, these state functions are the players 
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of multi-player game. A player's strategy must take into 
account of the other players's strategy. So game theory 
is born to be an entangled theory. This may help us to 
understand the entanglement among coexistence phases. 
We find a new quantity to measure the entanglement 
of different phases. The maximal point of entanglement 
is the Nash equilibrium point at which different phases 
coexist. The maximal entanglement between phases is 
controlled by interaction parameter instead of time, so 
the quantum entanglement between different quantum 
phases are good candidates for quantum computation. 

Further more, if we take two particles as the two play- 
ers of prisoner dilemma, it would be found that there are 
two types of entangled pairs: the angel pair and devil 
pair, the angel pair love each other, their entanglement 
is a little bit weaker. While the devil players do not trust 
each other, they inclined to bounded together to prevent 
the other from betray, it is a strong entangled pair. This 
may help us to understand pairing mechanism in super- 
conductor. 
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APPENDIX A: EHRENFEST'S DEFINITION 
ABOUT ORDER OF PHASE TRANSITION 



We rewrite Ehrcnfcst's definition into more compact for- 
mulism. The pth-order quantum phase transition is char- 
acterized by a discontinuity in the pth derivative of dif- 
ference of free energy SF, 

ee-HF = o, a^dsz m -HF = o, (m = i, 2, ..., P - 1), 



71 

d^SF = 0. 



o^sf + 0, d™d?; m SF + 0, (m = 1, 2, ..., P ), 

dPJF^O. (A5) 

If the pth derivative of 5F becomes continuous, the phase 
transition jumps to the p + lth order. 



APPENDIX B: DIFFERENTIAL GEOMETRY OF 
FREE ENERGY 

We demonstrate the basic conception of differential ge- 
ometry on the free energy manifold expanded by two 
thermal variables temperature T and pressure P. 

In differential geometry, the open set SF(T, P) in a 
manifold may be mapped to an open set in three di- 
mensional Euclidean space through the homeomorphic 
mapping f$ F : SF(T,P) — > X 5i r, X = Xi (71, 72)1 + 
^2(71,72)^ + ^3(71, 7 2 )/c, At each point 7 in SF there 
exist a set of tangent vectors. The two basis of the tan- 
gent vector space on this manifold are 



ei 



d 
^7i 



d_ 

<9T' 



e2 = 



d 
d-f2 



d_ 

<9T' 



(Bl) 



The zeroth order phase transition defines that the free 
energy of the two phases F(T, P) are not equal, 



F A (T,P)^F B (T,P). 



(Al) 



For the first order phase transition, the free energy of the 
two phases Fa and Fb is continuous, but the first order 
derivative is not continuous, 

TP IT P\- TP IT P\ dFA U dFB dFA U dFB 

F A (T,P)-F B (T,P), -^r^-^r, + 

(A2) 

The second order phase transition is defined by the dis- 
continuity of the compressibility, susceptibility, 

Cf^Cf, a A ^a B , k a ^k b , (A3) 

in which the specific heat C p , a,n arc defined as: 

r -r( dS ) r d2p 

_}_,9V_, 1 d 2 F 
a ~ V^df' P ~ V dTdP' 
l,dV. ld 2 F 

K = -v { dP )p = -vdP^- (A4) 

The higher order of phase transition is defined from the 
discontinuity of higher order derivative of free energy. 



here we have defined 71 = T(temperature) and 71 = 
P (pressure). An arbitrary tangent vectors may be ex- 
pressed as 

ax 9 2 x 

X 7 , - — , X 7i7j . - (B2) 

All the tangent vectors to the surface SF at x denoted 
form a tangent vector space denoted by T X SF, it is par- 
allel to the tangent plane to SF at x. 

The metric tensor is given in a bilinear form by the 
inner product of two basis vectors, gij 
the first fundamental form is 



< ej, ej >. Then 



E — ■ -^-71 ) F — -^-71 ' 


-^72 1 


G = 


5ii = E, 912 = 921 = F', 


922 = 


G. 


The second fundamental form is 






L — — 'X.™^™^ ■ 11 


-X 7l ■ 


n 7l 


M — X 7l72 • n = 


— X 7l 


' n 72 


N — X 7272 • n = 


— X 72 


• n 72 . 



'•72 -I 
(B3) 



with n as the normal vector of the surface, 



|X 71 X X 71 1 



(B4) 



(B5) 
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Considering the second order of phase transition, the first 
fundamental form and the second fundamental form in- 
cludes more information about the variation of the free 
energy than Ehrcnfest's definition. It is more reasonable 
to choose the Gaussian curvature, 



LN - M 2 
EG - F' 2 : 



(B6) 



It is positive for spheres, negative for one sheet hyper- 
bolic surface and zero for planes. The principal curva- 
tures, ki(7°) and ^2(7°), of X at X(7°) are defined as 
the maximum and the minimum normal curvatures at 
X(7°), respectively. The directions of the tangents of 
the two curves that are the result of the intersection of 
the surface X(7°) and the planes containing 11(7°). Then 

^(71,72) = ^rj(7i,7 2 )| p = ^]ki(7i, 72)^2(71, 72)| P - 
p p 

(B7) 

The cross section of the curvature X at 7 may be ex- 
panded as parabola by ignoring the higher order, 



X 



K l 2 . K 2 2 



(B8) 



where n = {dr/i + ^} J d r yid r yj)^/g^ + ■■ ■ . When k% > 
and k 2 > 0, it is a elliptic surface, when k\ > and 
k% = 0, it is parabolic, for n\ > and K2 < 0, it is 
hyperbolic. 

The covariant derivative may be defined as 



(B9) 



here we have introduced the Christoffel-Levi-Civita con- 
nection r^. , which is defined as 



ij 



1 JW/%2 , d 9. 



(BIO) 



APPENDIX C: DEFINITION OF GAME THEORY 

A game is an abstract formulation of an interactive 
decision situation with possibly conflicting interests, it 
involves of a number of agents or players, who are al- 
lowed a certain set of moves or actions. The payoff func- 
tion specifies how the players will be rewarded after they 
have performed their actions. Let n = 1, . . . , N denote 
the players; The ith player's strategy, Si, is her proce- 
dure for deciding which action to play, depending on 
her information. The strategy space of the zth player 
Sj={e,i, ei2, ■ ■ ■ , eim} is the set of strategies available to 
her. A strategy profile (si, S2, . . . , sat) is an assignment 
of one strategy to each player. Ui{s\, . . . , sn) is player 
i's payoff function (utility function), i.e., the measure of 
her satisfaction if the strategy profile (si, . . . , sn) gets 
realized. 

The strategy assigned to each player could be pure 
state or mixed state. A mixed strategy is a probabil- 
ity distribution over pure strategies. Playing a mixed 



strategy means, the players come up with one of her fea- 
sible actions with a pre-assigned probability. Each mixed 
strategy corresponds to a point P of the mixed strategy 
simplex. 

m 

M p = {P = (p!, . . . , Pm ) £ R m : Pq > 0, Y^Pi = 1}. 

9=1 

(CI) 

whose corners are the pure strategies. For the case of two 
players, a strategy profile is a pair of probability vector 
(p, q) with p 6 M p and q € M q . The expected payoffs 
of player 1 and 2 are expressed as 

m 

Wi(Pi q) = P • Aq = ^J/>,.1,. ,'/.;• (C2) 

ij 
m 

ua(p, q) = P Bq = ^2 Pl B ijqj . (C3) 

This payoffs function relies on the entangled strategies of 
the two players. 

The essence of a game is the payoff function defined 
over different strategy profile of the players. In the lan- 
guage of quantum statistics, the strategy space of the ith 
player Sa is the Hilbert space of the i particle. 

APPENDIX D: DIFFERENTIAL GAME 

The differential game is applied to model dynamic 
conflicts [HI, such as the labor-employers relations in eco- 
nomic processes, the bull and fighter in bull-fighting. The 
labor and employer are called players in differential game. 
We take the two player game as an example. There are 
two players, one is called Alice, the other player is Bob. 
Alice has a strategy space Sa(si, sa, • • • , sn), and Bob's 
strategy space is Sb- We introduce the state vector x 
which characterize the conflict to convert a real life con- 
flict to a differential game model. For example, Alice and 
Bob are concerned with the motion of a point in a plane, 
the state vector denotes the location of the point x(t) 
at the time t. All possible state vector is a subset of a 
rt-dimensional Euclidean space. 

The influence of the decisions of the evolution of the 
state is described by the equations of motion, 



/(x, s Al s B ) 



(Dl) 



When the two players realize their strategy (sa,&b) at 
time t, the outcome of the differential game is functional 
on the state space, J7[x, Sa, Sb]. If the outcome satisfies 
for a pair of strategy (s^,s B ) satisfies, 



<U[x,s A ,s* B ], (D2) 

then this strategy pair is the optimal play for the two 
players, U[x, s* A ,s* B ] is called the optimal outcome at 
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x. The optimal outcome is a function in the Euclidean 
space, it is called the value function, 



J(x) = U[x,s* A ,s* B ). 



(D3) 



There is a theorem states that if the value function of 
a differential game exists it is unique (38|. The saddle 
relation implies that a solution of a differential game is 
a Nash equilibrium. Neither player Alice nor player Bob 
can improve their guaranteed results. 



av(T,j>*) 



For 



The conjugate momentum of ib is p(T) 
a stochastic differential game, the Hamilton function is 



H s (t, i>, 7) = V^V(t, ^)/(t, 1>, 7) + g(t, V, 7), (E6) 



the Hamilton-Jacobi-Bellman equation reads 



dV(t,j>) , 1 y_^V" 

-4- — /T J 



+ minH s (t,ip,j) = 0. (E7) 



APPENDIX E: THE CONTINUOUS-TIME 
INFINITE DYNAMIC GAME 

We present the basic conception of the continuous-time 
infinite dynamic game in this section. Suppose for a given 
many-body system, there are n player which are denoted 
as N = {I, 2, n}, ib is the state function of this game, 
the state space is an entangled space of many Hilbert 
space. The evolution equation of this game is governed 
by the differential equation, 

*'^=/(*,V',7i(t),72(t),...,7n(«)). (El) 

71,72, ...,7n is the strategy profile of the n players. 74 are 
real or complex numbers. A cost function of the game is 
a map from the strategy space r = T\ ® T\ ® ... ® T n to 
a real number, E % : T — > R(i € AT) for each fixed initial 
strategy profile 7°, 7°, 7„- The minimum cost-to-go 
from any initial state and any initial time is described by 
the so-called value function which is defined by 



V(t, ib) — min 

7(*) 



g(s,m,l(t))ds + q(T^(T))}, (E2) 



satisfying the boundary condition V(T,ip) = q(T,ip), g 
is a map from g : (t,ip(t),^/(t)) — > R. The application of 
the principal of optimality leads to the Hamilton-Jacobi- 
Bellman equation, 



dV(t,iP) 
dt 



mm I 

7 



t d V(t,i/>) 
dip 



/(t,0,7)+5(t,^7)]- (E3) 



A theorem states that, if a continuously differentiable 
function V(t, ip) can be found that satisfies the Hamilton- 
Jacobi-Bellman equation, with corresponding boundary 
condition V(T,ip) — q(T,ip), then it generates the op- 
timal strategy through the static minimization problem 
defined by the right hand side of Eq. (|E3p . We introduce 
the Hamilton function 

H(t, V>, 7) = ^ 7) + g(Jt, 0, 7), (E4) 

the minimizing 7 will be denoted by 7*, then 

H(^, 7 *) + ^|^ = 0, (E5) 



For a given two person game, suppose that the strat- 
egy pair (71,72) provides a saddle-point solution, and 
ib* (t) denoting the corresponding state trajectory, then 
the Hamilton function satisfies, 



n(t, r , 71* , 72) < n(t, r , 71* , 72 ) < n(t, r , 71 , 72 * ), 

(E8) 

(71,72) is the solution of the fssac equation min max 7i = 
—dV/dt, V(tp) is the value function. It represents the 
coexistence hyper-surface extended in ib space. 



APPENDIX F: BRIEF INTRODUCTION TO 
COOPERATIVE GAME 



Coalition is a crucial ingredient in a n-person cooper- 
ative game. When a subset g of the n player forms a 
coalition and all its members act together, we assign a 
real number E(g) to each possible coalition, E(g) mea- 
sures the payoff when this coalition acts together. If we 
can find an imputation vector £ — {e\ 1 £2 1 ...,^n ) with 
real components, which satisfies £j > E({i}),\/i € and 
£1 + £2 + ■■■ + £ n = E(N), it is the Pareto optimality. 
An imputation x represents a realizable way that the n 
player can distribute the total payoff E(N). 

Usually, two coalitions S and T are inclined to united, 
if their union brings them better payoffs, this is the super- 
additivity of a game, E(SUT) > E(S) + E(T). Suppose 
a coalition S has m players, usually there is a difference 
between the payoff of the coalition and the sum of payoff 
of each individual in this coalition, the excess is defined 
by, 



e{S,e a ) = E{S)-J2 £ i 



(El) 



A n-person cooperative game has 2™ subsets, i.e., it has 
2™ coalitions. Let Z?2"-2 : M 2 "~ 2 — > M 2 "~ 2 be a mapping 
which arranges elements of a 2™ — 2-dimensional vector 
in order of decreasing magnitude. For a certain set e of 
payoff vector, the nucleolus over e is the solution mini- 
mizing the vector of the excesses, 
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N C (N, E, e) = {e a e e\D 2n _ 2 {e{S 1 ,e a ), e(5 2 »_ 2) e„))<£> 2 »-2(e(Si, e 6 ), e(S 2 »- 2) e 6 ))}, V6 G n. (F2) 

I 



In a game (JV, E) , the e-core C e (N, E) is the set of all 
pre-imputations e, satisfying that all the excess function 
are not greater than e, i.e., 

C e (JV,£) = {f a e e(N,E)\e(S,e a ) < e,VS C JV}. (F3) 

where the set of pre-imputations e(N,E) is e(N,E) = 
{e a eR n \j: i e ai = E(N),e ai >0}. 

The most important quantity we want to find is 
the ground state payoff vectors e g , i.e., the least 
core C e (N,E), which satisfies e(S,e g ) < e g = 
minmaxe(5, e). In fact, this is equivalent to the self- 
content mean field theory. The ground state payoff vector 
can be obtained by solving the following linear program- 
ming problem, which minimizes the maximal excesses: 



(F4) 

bject to E(S) - J2e l < e g ,VS C A, (F5) 



minimize e 



ei +e 2 + +s n = E(N), 
e i > 0,i = 1, 2, n. 



(F6) 
(F7) 



There is a lemma says the unique payoff vector e* mini- 
mizing s can always be determined by at most n steps in 
algorithm computation]!!]. 



APPENDIX G: PRISONER DILEMMA 

The Prisoner's dilemma is the most famous paradigm 
to study game theory. It occurs between two prisoners: 
Alice and Bob, who are accomplices to a crime which 
leads to their imprisonment. Each has to choose between 
the strategies of confession or accusation. If neither con- 
fesses, moderate sentences (a years in prison) are handed 
out. If Alice confesses and Bob accuses him, Bob is free 
(0 years in prison) and Alice is sentenced to c > a years 
in prison. If both confess, they will each have to serve b 
years in prison, where a < b < c. 







Bob 


Bob 






Not Confesses 


Confesses 


Alice 


Not Confesses 


(a, a) 


(c, 0) 


Alice 


Confesses 


(0, c) 


(b, b) 



If there is no cooperation and communication between 
the players, each of the players choose the strategies to 
minimize his losses as far as possible. We suppose the 
players Alice and Bob choose their strategies using their 
loss function ua and u B from Sa x Sb to R. The two 
players choose the strategies to minimize the biloss map- 
ping 



A consistent pair of strategies following the constrain 

/a(sb) := {sa G Sa\ua(sa, 8b) = inf u a (sa, s b )}. 
Ib{sa) ■= {sb € S b \u b (sb,s a ) = inf u b (sa, sb)}(G2) 

is called a non-equilibrium (Nash equilibrium). In the 
prisoner's dilemma, the Nash equilibrium is that both Al- 
ice and Bob confesses, they both spend b years in prison. 

Besides the Nash equilibrium, there is still a better 
strategies for them. If Alice and Bob communicate and 
cooperate with each other, they would not confess so that 
they only serve a < b years in prison. 



APPENDIX H: COURNOT DUOPOLY 

A market comprised of two sellers and many compet- 
itive buyers is know as a duopoly. The buyers can not 
influence the price or quantities offered, it is assumed that 
the collective behavior of the buyers is fixed and known. 
The competitive and cooperative behavior between the 
sellers determines the price. 

Two players are each manufacturers of the same single 
commodity, the loss functions are cost function which 
depends on the production of the two players. We denote 
the quantities of this commodity produced by the two 
players by x G K + and y G K + . The price p(x,y) is an 
affine function of the total production x + y, 



p(x + y):=a- 0(x + y), 



(HI) 



and the cost function ua and Ms of each manufacturer are 
affine functions of the production Jax = ax + b, /s(y) = 
ay + b. Alice's net cost is equal to /a(x, y) '■— j "a — p(x + 
y)x, Ja{x, y) := f B — p(x + y)y . Eliminating the constant 
terms which does not modify the game, the loss function 
reduced to 



u A (x,y) :=x(x + y-u), 
u B {x,y) := y(x + y-u). 



(H2) 
(H3) 



u(s A ,s B ) := {u a (sa,s b ),u b (s a ,s b )) G 



(Gl) 



The non-cooperative equilibrium may be attained algo- 
rithmically following the scenario of two payer's game. 



APPENDIX I: THE ATIYAH-SINGER INDEX 
THEOREM 

The Atiyah-Singer index theorem is concerned with the 
existence and uniqueness of solutions to linear partial dif- 
ferential equations of elliptic type. The Fredholm index 
is a topological invariant of elliptic equations. By com- 
puting a small number of fundamental examples and by 
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showing that both functions have similar algebraic prop- 
erties, Atiyah and Singer proved that the Fredholm index 
and the topological index are both topological invariants 
of elliptic equations, they are equal. 

According to the Atiyah-Singer index theorem [34l l35j 
the analytic index of the operator D is defined as [39] 



Index D = dimKerD — dimCokerD, 



(II) 



where Ker D is the kernel of the operator, which is defined 
to be the space of zero- modular solutions, i. e., on the 
entire space V 



KerD = {£ e T(E)\D£ = 0}. 



(12) 



On a real oriented compact smooth n = 2/-manifold M , 
the Atiyah-Singer index theorem states that the Euler's 
characteristic Ch(M) is the sum of the Betti number, 

Ch{M) = Index(A,d) = ^2(-l) p dim R H% R (T(M), R). 

p 

(13) 

It is just the Gauss-Bonnet-Chern theorem, on the four 
dimensional manifold, there are more interesting case. 
Such as the Dirac operator, there are k zero modes of 
fermions coupled to the k-instanton gauge field in the 
fundamental representation. Atiyah-Singer index theo- 
rem which states that the the index of the Dirac operator 
is minus the first Chern class, where v + is the number 
of positive chirality zero-modes and V— is the number of 
negative chirality zero-modes, i.e., 



IndexD — Ch = v + 



1 

2^ 



n. 



(14) 



M 



phase to Chern character. We define the 1-form of the 
eigenf unction as 



dMR(t))) = d Rk \ip(R(t)))dR k , 



(33) 



the wedge product of two one form corresponds to the 
conventional Berry phase, i.e., the first Chern character, 

B l =d{^(R(t))\Ad\i>(R(t))) 
= e ijk {d Rj iP(R(t))\d Rk ^(R(t)))d 2 R. (J4) 

Then the second order generalization is 

B 2 = d{ip\ A d\ip)Ad(ip\ A d\ip) 



ttjki (d Rt ipi | d Rj ipi ) (d Rk ip t | d Rl ip^cPR. 



(J5) 



It is natural to arrive at the pth order generalization, 

B p = d{i)\Ad\i/j)---Ad{il>\Ad\ip). (36) 



2p 



In mind of our new definition of the general Berry phase, 
the p-form expression of Eq. (|J6p is 



B l = II S ( E 3 ~ E i) d5E i A d6E * ( J7 ) 

Bf = J| 5(Ej - Ei)d5E % A dSEi A • • • d5E % A dSEi (38) 

According to Eq. (jJ2|) , the complete Chern character can 
be expressed into beautiful form, 



APPENDIX J: TOPOLOGICAL QUANTITY 

Chern character 

For a complex vector bundle on the base space-time man- 
ifold M, whose structure group is the general k dimen- 
sional complex linear group GL(k,c), there exists the 
Chern character, which is an invariant polynomial of the 
group GL(k, c), 



Ch(M) = Tr(exp-n) 

- k+ i Trn+ kj^ TrinAQ)+ --< 31) 

While this invariant polynomial may be expressed into a 
more familiar generalized Berry phase, 



Ch B (M) = B + B 1 + j [ B 2 



(32) 



My Generalization of Berry Phase 

In this section, we generalize the conventional Berry 



Ch(M) = Tr[exp(-S(E j - E z )dSE, t A dSEi)} 

- k+ i B ^h(^r^ + -- (J9) 

When the correction to the zth energy level SEi is ex- 
panded up to the pth order, in mind of Eq. (??), it is 
easy to verify that 

dSEAdSEA ■ ■ ■ AdSEAdSE 

m 

= d8E {m ^hd5E^ m) h ■ ■ ■ Ad5E^Ad8E {m \ (J10) 

This equation suggests that if we want to find out the 
pth order of phase transition, the highest order of per- 
turbation to the ith energy level must be extended to the 
pth order. 
Berry Phase 

For the time dependent Hamiltonian operator Hi(R(t)), 
the eigenvalues and eigenvectors at time t can be ex- 
pressed as a function of t, Ei(R(t)) and tpi(R(t)). 

Hi(R(t))\ipi(R(t))) = E l (R(t))\ip i (R(t))} . (311) 
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According to Berry's definition [32|], there was a induced 
gage potential 



A(R(t)) = i(iPi(R(t))\VMi(R(t))). 



(J12) 



The first Chern number of Berry phase can be rewritten 
by these coordinates in R space as 

1 



Ch! = — I Vx A(R(t)), 

2lTT J R 

= ie ijk (d Rj xP(R(t))\d Rk 4>(R(t))). (J13) 

The Chern number is actually the topological quantiza- 
tion of the magnetic field 

B 1 = Vx A= ie t jk{d R:l i/Ji\d Rk ipi) 

{^\d Rt H\^ p ){^ p \d R] H\^) 



ie tJ k 



-Imetjk £ 



{E p - E,Y 



(J14) 



There has been some evidences pointing out that the 
Berry curvature is singular at the points where the energy 
bands touches [33(. 



APPENDIX K: LINEAR RESPONSE AND 
GREEN FUNCTION 

Let Hq be the full Hamiltonian describing the system 
in isolation. One way to test its properties is to couple the 
system to a weak external perturbation and to determine 
how the ground state and the excited states are affected 
by the perturbation. The Hamiltonian H for the system 
weakly coupled to an external perturbation, which we 
will represent by a Hamiltonian H e , is H = Hq + H e . 
Let 0(x, t) be a local observable, such as the local den- 
sity, the charge current, or the local magnetization. The 
expectation value of the observable in the exact ground 
state |0), under the action of the weak perturbation H e , 
is modified as 

(0|6|0) e = (0|O|0)+ l n- 1 / dt'(0\[H e ,6]\0)+---{Kl) 



If we only retain the linear term in H e , the first-order 
change in a matrix element arising from the external per- 
turbation is expressed in terms of Heisenberg operator of 
the interacting, S(0\6\0) = (0|6|0) e - (0|6|0), i.e. 



5(0\6\0) = ih- 1 dt'(0\[H e ,6]\0). 



(K2) 



This change represents the linear response of the system 
to the external perturbation. It is given in terms of the 
ground state expectation value of the commutator of the 



perturbation and the observable. If O is a local observ- 
able, H e (t) represents an external source which couples 
linearly to the observable, H e = J Of(x,t). The coeffi- 
cient of proportionality between the change in the expec- 
tation value (0\O(x,t)\0) and the force f(x,t) defines a 
generalized susceptibility 



dTe iuJT (0\[6(x',t'),6(x,t)}\Q). (K3) 



G(x' , x) = [0(x' , t'), 0(x, t)] is the retarded Green func- 
tion. Under Fourier transforms, this Green function can 
be mapped into its momentum space, 

G(k, u) = I d(x - x 1 ) I dte tk(x '- x) e lut G{x', x) a (K4) 



In fact, the correction to the ground state energy can 
also be expressed by Green function in momentum space. 
Usually, it is easy to get SE g from the familiar formula, 



5E n 



\He\0) 



|g e |rc)Mff e |0) 
Eq — E n 



(K5) 



We can write the Hamiltonian with a variable coupling 
constant A as H(X) = Hq + XH e ,so that H(l) — H and 
H(0) — Ho- Then the correction to the ground state 
energy is, 



V 



1 dX 



SE = ±-i- 

2 (2tt4) J q X 



d A k 



due^ij 



2 J.2 



2m 



)TrG (k, ui). 



(K6) 



Therefore, the topological quantum phase transition is 
intimately related to the topology of momentum space. 
In fact, the correction to the ground state energy is not 
the only physical quantity that can be used to study the 
quantum phase transition, other physical observables also 
present a sudden change at the phase transition point. 
The thermal dynamic potential is another good candi- 
date, 



5CI = fi - fi = ± 



1 dX 



[ d 3 zlimlimi[-^ + — + fj,]TrG x (xT, x'tftZ) 
J 2 or 2m 

There also exist many other physical observables, such 
as the single-particle current operator, J = J d 3 xj, 
Ja/3 — J^ap^aJapipfi, whose expression in terms of 
Green function is < j >— ±i \im\imTr[J(x)G(xt, x'tf)]. 
The number density operator of particles is < p(x) >= 
±iTrpG(xt, x't'). The spin density operator < a(x) >= 
±iTraG(xt,x't'). 
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APPENDIX L: GREEN FUNCTION 

We first show that the Green function follows a sim- 
ilar Schrodinger equation. For a system described by 
a Hamiltonian H which can not be solved exactly, the 
usual approach is to set, H = Hq + V, where Ho is the 
unperturbed part which may be solved exactly. The term 
V represents all the interactions. The wave function is 
governed by the interaction, 



(LI) 



The wave functions of the Schrodinger equation, 
idtip = Hip, are time dependent, ^(t) = U(t)ip(0) = 
e lHot e~ iHt ip(0).The operator U{t) obeys a differential 
equation which can written in the interaction represen- 
tation, id t U{t) = V(t)U(t). The operator U(t) in the 
interaction representation has a expansion, U(t) = 1 + 

En=l ^ Jo ^ /o.*2- ' '*» /o nV{tl)V{h) ■ ■ ■ V(tn)}, 

T is the time-ordering operator. U (t) may be abbre- 
viated by writing it as U(t) = Texp[— i J Q dtiV(ti)].. 
The S matrix changes the wave function ip(t') into tp(t) 
may be defined from U(t), $(t) = S(t,t')ip(t'). The 
time-ordered operator also obeys S(t,t') — U(t)U(t'). ft 
is easy to verify S(t,t') obeys 



id t S(t,t') = V(t)S(t,t'). 



(L2) 



This equation is intimately related to the computation 
of ground state (or vacuum) expectation values of time 
ordered products of field operators in the Heisenberg rep- 
resentation 



The terms in the series for (0|5(+oo, — oo)|0) are called 
vacuum polarization terms. There was theorem states 
that the vacuum polarization diagrams exactly cancel the 
disconnected diagrams in the expansions for G(p, t — t'). 

A great simplification is that we sum over only 
linked (connected) distinct graphs. It turn out that 
this corresponds exactly to cancelling the factor of 
(0|S , (+oo, — oo)|0) in the denominator. From the discus- 
sions above, one can that the unperturbed part of the 
Hamiltonian Hq plays a trivial role in the quantum phase 
transition, it is the interaction part H e that decide the 
orientation of the evolution of a physical system. If we 
introduce the self-energy function T,(p, E) to absorb all 
the interaction, the exact Green function can obtained 
from the Dyson's equation, 

Gip^y 1 ^G^(P,E)[l-G ( P ,m(p,E)}. (L7) 

we expand G in the power series G(p, E) = Gq + Gq 2 ^ + 
Go 3 S 2 + • • • . The self-energy is a summation of an infi- 
nite number of distinct diagrams in the series. However it 
is impossible to get E(p, E) exactly, one must be contend 
with an approximation results. Usually, the higher order 
of phase transition means the higher order of self-energy 
terms must be included. The Green function of an inter- 
acting electron system is G r (p, uj) = [u> — Ek — S] ~ 1 . The 
corresponding spectra function is A(p, u>) = ——ImG r , 



A(p,w) 



lim 

7mS r - 



T i 'wImS r 



o ( w - E'Y 



5(uj-E'). 



G N (x u x 2 ,...,x N ) = (0\T[<j>( Xl )cj)(x 2 )...<l>(x N )]\Q), (L3) where E' g = E q - ReS r (p,u). 



(L8) 

The spectra function 



In particular the 2-point Green function, Gr{x\,X2) = 
(0\T{(f)(xi)(f)(x2)]\0) , is known as the Feynman Propaga- 
tor. For an interacting case, the Green function reads [39j| 

12 _ (0|T[^ 1 )^ 2 )g(+oo,-(x j )]|0) 

U (Xi,X 2 — , . , rr-r • (1j4J 

(0|5(+oo, — oo)|0) 

The Green function of energy is defined by the usual 
Fourier transformation with respect to the time invari- 
able: 



represents a resonance peak with width 2r 9 . One can 
associate each peak with a quasi-particle. The life- 
time of quasiparticle is infinite for a vanished JmS r , 
T q = lim/ m s r [2JmE r ] _1 — > oo. So the quasiparticles on 
the fermi surface is a kind of topological excitation, the 
external perturbation does not shorten its life. As all 
know, E = Gq(P, E)^ 1 is the unperturbed part of the 
system. Therefore it is the self-interacting part S(p, E) 
that determines a phase transition, in the frame work of 
out theory, one can denote the self-energy as the pertur- 
bation part to exactly solved part, 



/-t-oo 
dte iE ^G(p,t-t'). (L5) 8E> = G(p,E)- l -G Q {P,E)- l = V{p,E). 

-OO 



(L9) 



According to the Goldstone's theorem, 

0|ff e S'(+oo,-oo)|0) 



SE — E — Eq 



(0|S(+oo,-oo)|0) 



(L6) 



Usually E(p, E) is a complex matrix, it may be rewrit- 
ten in the Dirac's bra and ket representation, S(p, E) = 
\X(p,E)) and £(p,£)t = 
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